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PREFACE 


This  report  was  prepared  by  the  University  of  Dayton 
Research  Institute  for  the  Federal  Aviation  Administration 
Systems  Research  and  Development  Service  under  Contract  FA74WA- 
3532  during  the  period  June  1976  to  December  1977.  The  report 
describes  a comparison  of  Dayton  Aircraft  Cabin  Fire  Model 
(DACFIR)  to  a series  of  aircraft  cabin  mock-up  fire  tests  and 
the  refinement  and  extension  of  the  model  based  upon  the 
comparison. 

Technical  administration  of  this  work  was  provided  by 
Mr.  Charles  C.  Troha  of  the  Systems  Research  and  Development 
Service  (ARD-520) . Work  was  performed  at  the  University  of 
Dayton  under  the  supervision  of  Mr.  Nicholas  A.  Engler. 

Other  personnel  at  the  University  who  have  contributed  to  this 
program  include  Mr.  Peter  M.  Kahut,  Mr.  Steven  G.  Vondrell, 
and  Mr.  Richard  E.  Feldmann.  The  laboratory  data  collection 
program  performed  under  subcontract  by  the  Boeing  Commercial 
Airplane  Company  was  directed  by  Dr.  James  M.  Peterson  with 
contributions  by  Mr.  Robert  D.  Lechner,  Dr.  Allen  E.  Senear, 
and  Mr.  Everett  A.  Tustin.  Much  of  the  original  development 
of  the  DACFIR  model  was  directed  by  Mr.  Jerry  B.  Reeves  who 
also  guided  the  early  work  on  this  refinement  program.  The 
authors  would  like  to  thank  Ms.  Jacquelin  Aldrich  for  her  patient 
assistance  in  preparing  the  manuscript. 
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SECTION  1 
INTRODUCTION 

This  report  describes  the  results  of  an  evaluation  and 
refinement  of  the  Dayton  Aircraft  Cabin  Fire  (DACFIR)  Model. 

The  DACFIR  model  was  used  to  make  computer  simulations  of 
several  full  scale  aircraft  cabin  and  cabin  mock-up  fire  tests. 
Comparisons  of  the  test  results  to  the  model's  predictions  were 
used  to  verify  the  computer  code,  evaluate  the  model's  performance, 
and  refine  the  model. 

The  DACFIR  model  was  developed  by  the  University  of  Dayton 
Research  Institute  (UDRI)  for  the  Federal  Aviation  Administration 
(FAA) , to  enable  the  smoke  and  toxic  gas  emissions  from  the 
burning  of  interior  materials  in  a cabin  fire  to  be  predicted 
from  laboratory  test  data  on  these  materials.  The  development 
of  the  basic  mathematical  model  is  described  in  Volumes  I,  II, 
and  III  of  Department  of  Transportation  (DOT)  Report 
FAA-RD-76-120  . 

Originally  the  DACFIR  model  was  developed  specifically  for 
simulating  fires  in  wide-body  aircraft.  Certain  parts  of  the 
original  computer  program  concerned  with  the  cabin  geometry 
limited  the  model's  application  to  this  type  of  aircraft.  No 
full  scale  fire  tests  have,  however,  been  conducted  for  wide- 
body  cabins  while  the  results  of  several  series  of  tests  of 
standard  body  cabins  are  available.  In  order  to  evaluate  the 
model's  performance,  a program  was  initiated  to  do  the  following. 

(1)  Modify  the  DACFIR  program  so  that  both 
wide-body  and  standard-width  cabin 
geometry  could  be  used. 

(2)  Collect  information  on  past  full  scale  burn 
tests  and  select  a number  of  these  tests  to 
be  simulated  with  the  DACFIR  model. 


^ Reeves,  J.B.  and  C.D.  MacArthur,  "Dayton  Aircraft  Cabin  Fire 
Model",  Volumes  I,  II,  and  III,  FAA-RD-76-120,  June  1976. 


1 


(3)  Collect  laboratory  test  data  on  the 
interior  materials  used  in  the  selected 
full  scale  tests. 

(4)  Simulate  each  full  scale  test  with  the 
DACFIR  model  and  compare  the  model's 
results  to  those  of  the  test. 

(5)  Identify  areas  in  which  the  DACFIR  model 
requires  refinement  and  make  the 
indicated  refinements  where  possible. 

Based  upon  the  results  of  the  comparison  and  refinement 
process,  suggested  guidelines  were  established  for  future  full- 
scale  tests  to  validate  the  model. 

Modifications  made  to  the  program,  both  for  the  purpose  of 
simulating  fires  in  standard  body  cabins  and  as  a result  of  com- 
parisons to  full  scale  tests,  are  described  in  Section  2 of  this 
report.  Section  3 presents  the  results  of  the  comparison  of  the 
refined  model  to  the  full-scale  tests  and  analyzes  the  model's 
performance  in  simulating  these  cases.  Section  4,  Conclusions, 
summarizes  the  results  of  this  study.  Four  appendices  are  added 
to  the  report.  Appendix  A presents  the  development  of  a proposed 
major  refinement  to  the  DACFIR  model  concerning  the  simulation 
of  the  cabin  atmosphere  dynamics.  Appendix  B describes  the 
laboratory  data  collection  program  conducted  in  support  of  the 
validation  exercise.  Appendix  C is  a User's  Guide  for  the  computer 
program  which  implements  the  refined  version  of  the  model. 

Appendix  D contains  derivations  of  several  equations  presented 
in  the  body  of  the  report. 
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SECTION  2 

REFINEMENT  OF  THE  DACFIR  MODEL 

This  section  describes  the  modifications  and  refinements 
of  the  DACFIR  model  made  as  a result  of  this  validation  exercise. 

The  modified  and  refined  version  of  the  model  is  designated  as 
Version  2 of  DACFIR  or  DACFIR2  for  short.  DACFIR2  retains  the  basic 
mechanism  for  the  representation  of  the  fire  involvement  of  the 
cabin  interior  materials  used  in  the  original  model.  The 
paragraphs  below  give  a brief  review  of  the  structure  of 
the  DACFIR  model.  For  a more  complete  description,  see 
Reference  [1] . 

In  the  DACFIR  model,  the  interior  surfaces  of  an  aircraft 
cabin  - the  floor,  ceiling,  sidewalls,  etc.  - are  assumed  to  be 
flat  and  have  either  an  exactly  horizontal  or  vertical  orientation. 
To  chart  the  progress  of  a fire  burning  on  or  impinging  on  a 
surface,  the  surface  is  divided  into  square  elements  0.5  feet  on 
a side;  the  fire  behavior  of  the  material  composing  the  surface 
being  approximated  by  monitoring  the  state  of  each  element.  Four 
"fire  behavior"  states  are  assumed:  (1)  virgin;  (2)  smoldering; 

(3)  flaming;  and  (4)  charred.  An  element  of  material  in  the 
virgin  or  charred  state  does  not  emit  smoke,  heat,  or  gases 
while  an  element  in  the  smoldering  state  may  emit  smoke  and 
gases  but  not  heat.  Elements  in  the  flaming  state  are  the  active 
participants  in  the  fire  emitting  heat,  smoke,  and  gases.  Further, 
the  rates  of  heat,  smoke,  and  gas  release  for  flaming  elements  are 
functions  of  the  imposed  heat  flux,  mainly  radiant,  fed  back  from 
the  flames  of  the  fire  involving  the  element. 

Transitions  of  elements  from  one  state  to  another  occur  by 
several  mechanisms:  creeping  flame  spread  over  a surface  from 
groups  of  flaming  elements  to  adjacent  elements  not  yet  ignited, 
contact  and  envelopment  by  flames  from  a near-by  fire  striking 
a surface,  and  the  transition  to  smoldering  caused  by  the  radiant 
level  from  a near-by  fire.  The  rates  and  times  that  govern  these 
transitions  as  well  as  the  rates  of  emission  of  heat,  smoke,  and 
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gases  are  quantities  supplied  as  input  data  for  the  program  and 
are  obtained  from  laboratory  measurements  made  on  small  samples 
of  the  cabin  materials. 

Once  the  DACFIR  model  has  determined  at  a given  time  in 
the  simulation  the  states  of  all  elements  and  the  rates  of 
emission  of  heat,  smoke,  and  toxic  gases  from  those  elements 
actively  or  passively  involved  in  the  fire,  the  next  task  is 
to  determine  the  condition  of  the  cabin  atmosphere.  The  emitted 
heat,  smoke,  and  gases  are  used  as  input  to  a modal  of  the  cabin 
atmosphere  to  update  its  descriptive  parameters:  temperature, 
visibility,  and  gas  concentrations  which  bear  upon  occupant 
survival  and  escape.  The  cabin  gas  dynamics  model  divides 
the  cabin  atmosphere  into  two  zones,  an  upper  zone  consisting 
of  combustion  products  and  heated  and  vitiated  air,  and  a lower 
zone  consisting  of  cooler  ambient  air.  The  fire  pumps  lower 
zone  air  into  the  upper  zone  which  grows  in  thickness  at  a rate 
determined  by  the  fire  size,  the  loss  of  gas  through  cabin 
exits,  and  the  transfer  of  heat  from  this  zone. 

The  following  sections  present  the  refinements  made  to  the 
basic  model.  To  simulate  fires  in  cabins  of  standard  width, 
added  flexibility  was  provided  in  the  model's  description  of 
cabin  geometry.  The  gas  dynamics  calculations  were  upgraded 
so  that  oxygen  depletion,  forced  ventilation,  and  the  effect 
of  a circular  cabin  cross-section  could  be  included.  The 
radiation  heat  transfer  computations  were  refined  by  adding  a 
gray  gas  approximation  for  upper  zone  radiation  and  by  relating 
the  flame  radiation  to  the  flame  soot  concentration. 

An  alternative  approach  to  modeling  the  cabin  atmosphere, 
in  which  calculation  of  fore-aft  temperature  and  concentration 
gradients  are  made,  was  investigated  as  a candidate  replacement 
for  the  present,  spatially  lumped  gas  dynamic  model.  While  this 
new  approach  appears  promising  it  has  not  as  yet  been  sufficiently 
developed  so  that  it  may  be  integrated  into  the  complete  DACFIR2 
model.  For  this  reason,  a description  of  this  horizontal- 
gradient  model  has  been  placed  in  Appendix  A of  this  report. 
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2.1  CABIN  INTERIOR  GEOMETRY 

Version  1 of  DACFIR  was  constructed  specifically  for  the 
simulation  of  fires  within  wide-body  aircraft  cabins.  The  wide- 
body  cabin  arrangement  of  seats,  overhead  compartments,  and 
other  features  does  not  vary  significantly  between  different 
wide-body  models.  Sufficient  flexibility  was  incorporated  into 
Version  1 so  that  all  of  the  current  wide-body  cabin  configur- 
ations could  be  handled.  However,  no  full  scale  burn  tests  of 
wide-body  cabins  were  available  while  a relatively  comprehensive 
and  well  documented  amount  of  test  data  could  be  found  for 
cabins  of  standard  width,  i.e.,  that  of  the  B 707/727/737  or 
DC-8/9  width.  Modifications  were  therefore  required  to  allow 
the  DACFIR  model  to  represent  cabin  configurations  for  standard 
body  widths  including  the  older  overhead  structure  of  shelf- 
like hatracks  and  the  more  curved  sidewall  and  ceiling  surfaces 
of  these  smaller  fuselages. 

A technique  was  developed  through  which  the  user  of  the 
DACFIR2  program  need  only  supply  a few  overall  dimensions  of 
the  cabin  lining  surfaces  (i.e.,  the  floor,  sidewalls,  hatracks/ 
stowbins,  etc.)  and  the  components  of  the  unit  vector  normal 
to  the  surface.  The  program  assembles  the  surfaces  in  the 
proper  arrangement,  divides  the  surfaces  into  the  unit  elements, 
and  sets  all  indices  and  counters  necessary  for  the  fire 
simulation  to  proceed  over  the  surfaces.  All  common  transport 
category  cabin  interior  linings  can  be  represented  by  this 
method . 

The  increased  geometric  flexibility  required  a more  rigorous 
definition  of  the  cabin  coordinate  system.  All  dimensions, 
locations,  and  directions  in  DACFIR2  are  now  specified  in  a 
single  coordinate  system.  A right-handed  cartesian  system  is 
used  with  the  origin  located  in  the  forward  lower  right-hand 
corner  of  the  cabin  as  viewed  facing  forward.  This  coordinate 
system  is  shown  in  Figure  2-1.  The  figure  also  shows  a typical 
location  for  the  "detailed  section",  that  part  of  the  cabin 
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Figure  2-1.  Cabin  Coordinate  System.  CL  is  the  cabin  length, 

Cw  is  the  cabin  width,  Cy  is  the  cabin  height,  Sl 
is  the  detailed  section  length.  The  detailed 
section  contains  the  materials  involved  in  the  fire 
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Figure  2-2.  Seat  Group  Configuration.  S^  is  the  seat  group 

width  given  as  a multiple  of  0.5  feet.  All  other 
seat  dimensions  are  fixed  as  shown  in  units  of  feet. 
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in  which  the  fire  is  assumed  to  originate  and  thus  where  the 
materials  are  divided  into  elements  for  the  tracking  of  the 
fire.  It  is  in  this  area  that  all  of  the  burning  is  assumed 
to  occur  during  the  simulation  period  and  therefore  a considerable 
savings  in  computer  memory  can  be  realized  by  keeping  track  of 
elements  only  within  this  region.  The  gas  dynamics  model, 
however,  considers  the  entire  cabin  volume. 

In  DACFIR2,  surfaces  are  assumed  to  be  planar  and  oriented 
horizontally  or  vertically  as  in  the  original  model.  The  vector 
normal  to  a surface  is  input  to  identify  the  surface  orientation 
and  a "z  displacement"  is  also  used  to  locate  the  bottom  edge 
of  the  surface  with  respect  to  the  cabin  floor.  For  horizontal 
surfaces,  the  z displacement  is  the  distance  from  the  floor  to 
any  point  on  the  surface  (as  all  points  are  equidistance  from 
the  floor) . A surface  width  is  specified  which  is  the  surface 
a dimension  for  vertical  surfaces  and  the  x dimension  for 
horizontal  surfaces.  The  dimension  of  all  surfaces  in  the  y 
direction  is  taken  to  be  the  same  as  the  detailed  section  length. 

Seat  positions  and  seat  row  widths  are  user  definable  in 
DACFIR2.  Seats  are  modeled  in  DACFIR2  with  the  same  bench-like 
approximation  of  the  original  model.  Seats  side  by  side  in  a 
row  are  regarded  as  a single  bench  seat,  called  a seat  group, 
whose  wi,dth  may  be  specified  provided  that  it  is  an  integer 
multiple  of  0.5  feet.  The  other  dimensions  of  a seat  are 
assumed  constant.  The  seat  cushion  is  1.5  feet  from  back  to 
front;  the  backrest  3.5  feet  high,  and  both  cushion  and  back- 
rest are  0.5  feet  thick.  Figure  2-2  shows  the  seat  configuration 
and  dimensions.  Up  to  nine  seat  groups  may  be  represented  in  the 
model,  each  with  a different  width  if  required.  All  seats  are 
assumed  to  face  forward,  the  vertical  backrest  planes  perpendicular 
to,  and  the  horizontal  cushion  planes  parallel  to  the  floor. 

The  location  of  each  seat  group  is  given  by  specifying  the  x and 
y coordinates  of  the  forward  right-hand  corner  of  the  seat 
cushion  as  viewed  by  a person  sitting  in  the  seat. 
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An  improved  method  of  specifying  cabin  doors,  window  exits, 
or  other  vents  in  the  cabin  walls  has  also  been  included  in 
DACFIR2.  In  the  first  version  of  the  program,  only  vents  whose 
opening  extended  from  floor  level  to  a specified  height  were 
allowed,  that  is,  only  exit  doors  were  recognized.  The  full-scale 
tests  selected  to  validate  the  model  used  both  natural  and 
forced  ventilation  from  air  ducts  at  various  heights  and  through 
doors  of  various  widths  and  heights.  To  accommodate  such 
features,  DACFIR2  requires  that  the  height  and  width  of  a vent 
be  given  and  the  distance  from  the  top  of  the  vent  to  the  floor  be 
specified. 

2.2  THERMAL  RADIATION 

Two  refinements  to  the  methods  of  computing  thermal  radiation 
heat  transfer  have  been  made.  The  first  concerns  radiation  from 
the  hot  gas  layer  to  the  surfaces  in  the  upper  and  lower  zones. 

The  second  is  an  improvement  of  the  estimation  of  the  local 
radiation  from  a fire  to  its  fuel  bed. 

In  the  original  DACFIR  model  the  upper  zone  gas  was  assumed 

to  be  opaque  and  radiate  as  a black  body  at  the  upper  zone  gas 

temperature  (see  Section  5.4,  pp.  53-56  of  [1]).  A step  away 

from  this  simple  estimate  of  the  radiation  has  been  taken  by 

(21 

adopting  the  method  of  Quintiere 1 J in  estimating  the  upper  zone 
gas  emissivity  as 

(2-1)  e ■ 1 - exp (- (kgS  + k )L] 

where  eg  is  the  emissivity  of  the  upper  zone  gas, 

k is  the  total  absorption  cross  section  of  the 

s soot  (smoke)  in  the  gas  (ft2/particle) , 

kg  is  the  gas  band  absorption  coefficient  (ft-1) , 

S is  the  upper  zone  smoke  concentration 
(particles/ft^) , and 

L is  the  upper  zone  thickness  (ft)  which  is  used 
as  the  radiation  mean  beam  length  in  this  case. 

f 2 1 

1 Quintiere,  J.,  "The  Growth  of  Fire  in  Building  Compartments", 
presented  at  the  ASTM-NBS  Symposium  on  Fire  Standards  and  Safety, 
National  Bureau  of  Standards,  Gaithersburg,  Maryland,  1976. 
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The  value  of  kg  was  chosen  as  0.1054  ft  /particle  of  smoke.  It 
is  assumed  that  the  soot  radiation  is  not  a function  of  wave- 
length (gray  gas)  so  that  the  absorption  by  the  soot  is  the  same 
for  thermal  radiation  as  it  is  in  the  visual.  Quintiere ' s value 
for  the  gas  band  absorption  coefficient,  k^,  of  0.1  ft  1 has 
been  adopted.  With  the  upper  zone  total  emissivity, e defined 
by  Equation  (2-1)  the  net  radiation  loss  terms  for  the  upper  zone 


gas  are 

(2-2a) 

qru 

= e oA  F (T  4 - T 4) 
g u u'  u su 

and 

( 2 —2b ) 

qrl 

— 4 4 

= e aA  F.(T  - T , ) 

g xs  1 u si  ' 

where 

qru 

is  the  net  radiation  exchange  between  the  upper 
zone  gas  and  the  upper  walls  and  ceiling  at 
temperature  Tgu, 

qrl 

is  the  net  radiation  exchange  between  the  upper 
zone  gas  and  the  lower  walls  and  floor  at 
temperature  Tg1 , 

A 

u 

is  the  upper  wall  and  ceiling  area. 

A 

xs 

is  the  area  of  the  interface  between  the  upper 
zone  gas  and  the  lower  zone  gas. 

F 

u 

is  the  effective  view  factor  for  the  exchange  , 

F1  is  the  effective  view  factor  for  the  gas  and 
lower  surfaces . 

Equations  (2-2a)  and(2-2b)  given  here  replace  Equations  (5-16c) 

r _ 

and  (5-17c)  of  [1] . The  expressions  used  for  Fu  and  F^  are  given 
in  Appendix  D. 

A second  refinement  of  radiation  heat  transfer  in  the  DACFIR 

model  concerns  the  computation  of  flame  radiation.  The  best 

estimate  of  the  radiation  feedback  from  a fire  to  its  fuel 

is  critical  in  the  DACFIR  model  since  this  radiation  1 •■vel  is 

used  to  select  the  proper  values  of  flame  spread  and  emission 

rates  from  the  material's  data.  The  original  model  used  formula- 

f 3 1 

tions  of  Dayan  and  Tien  for  the  radiation  from  a cylindrical 

flame  to  its  base  center  and  base  edge.  These  formulations 


^Dayan,  A.  f and  C.L.  Tien,  "Radiant  Heating  from  a Cylindrical  Fire 
Column",  Combustion  Science  and  Technology,  Vol.  9,  (1974),  pp. 41-47. 
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involved  the  assumpcion  of  a value  for  a , the  total  emittance 

G 

of  the  flame  at  the  center  of  its  base.  A simple  estimate  of 
this  quantity  was  made  as  a linear  function  of  the  fire  base 
radius  in  the  original  model.  In  DACFIR2  the  development  of 
Dayan  and  Tien  has  been  followed  further  to  relate  ac  not  only 
to  the  fire  base  radius  but  also  to  the  flame  absorption 
coefficient. 


For  sooty  flames,  the  type  expected  in  cabin  fires,  practically 
all  of  the  radiation  is  due  to  the  soot  particles  in  the  flame. 

Since  measurements  of  the  smoke  generated  by  the  cabin  materials 
are  available,  and  if  it  can  be  assumed  that  the  flame  is  a gray 
gas,  the  emissivities  of  fires  of  the  materials  can  be  computed 
from  the  smoke  generation  data.  Using  these  assumptions,  a 
method  was  developed  to  express  the  absorption  coefficient  in 
terms  of  the  smoke  generation  rate  as 


(2-3) 

where 


kf  = 0.21$" (hfg)-1/2 

kf  is  the  flame  absorption  coefficient  (ft  1) , 

is  the  smoke  generation  rate  of  the  fire 
fuel  (particles/ft2  . sec) , 

hf  is  the  flame  height  (ft) , and 

g is  gravity  (ft/sec^) . 


The  derivation  of  Equation  (2-3)  is  given  in  Appendix  D.  Using 
this  estimate  of  the  flame  absorption  coefficient,  the  emittance 
at  the  fire  base  center  is  [3] 


(2-4)  ac  = [1  - exp(-1.8kfhf ) ] 

+[1  - exp(-1.8kfyQ) ] 

-{I  - exp [-1 . 8kf (hf 2 + y02)1/2]> 

where  yQ  is  the  flame  cylinder  base  radius  (ft.).  Equation  (2-4) 
replaces  Equation  (4-2)  of  [1]  . 


Reconsideration  of  the  expression  used  in  the  original  model 
for  the  radiation  to  the  base  edge  led  to  a simplification 
by  assuming  the  appropriate  radiation  intensity  to  be  that  at 
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the  immedidate  edge  of  the  fire.  The  revised  expression  used 
in  DACFIR2  which  replaces  Equation  (4-1)  of  [1]  is 

(2-5)  q1  = 0.5ebac 

2 

where  is  the  assumed  flame  emissive  power  (16.3  Btu/ft  • sec). 
2.3  OXYGEN  DEPLETION 

The  depletion  of  oxygen  from  the  upper  zone  gas  by  the 

burning  of  the  cabin  materials  has  been  included  in  DACFIR2 . 

The  information  to  do  this  is  available  from  the  input  data 

items  for  each  material:  Q , the  heat  of  combustion  (Btu/lbm) , 

c 

and  y,  the  stoichiometric  oxygen  to  fuel  mass  ratio.  Using 
these  quantities,  the  heat  released  per  unit  mass  of  oxygen 
burned  is  given  by 

(2-6)  Fq  = Qc/Y 

where  F0  is  the  oxygen  consumption  factor  (Btu  released/ 

u lbm  02  consumed)  for  a given  material. 

Then,  using  the  heat  release  rate  for  the  material,  the  oxygen 
consumption  rate  is  given  by 

(2-7)  wQ  = q"/F0 

2 

where  wQ  has  units  of  lbm  of  02/ft  -sec  and  q"  is  the  heat 
release  rate  of  the  material  in  units  of  Btu/ft^-sec. 

At  each  time  step  in  the  simulation,  the  individual  oxygen 
consumption  rates  for  each  burning  material  are  computed  from 
Equation  (2-7)  using  the  appropriate  heat  release  rates  and 
areas  in  much  the  same  fashion  as  the  toxic  gas  generation  is 
computed.  The  gas  from  which  the  oxygen  is  removed  is  that 
forming  the  upward  flowing  fire  plume  so  that  the  upper  zone  is 
the  oxygen  depleted  region.  A conservation  equation  is  used  to 
compute  the  concentration  of  oxygen  in  the  upper  zone  gas, 

a?  <Y0Mu>  - ^Y0a*f  - V - Vo  * Vi 

where  YQ  is  the  mass  fraction  of  oxygen  in  the  upper  zone, 

Mu  is  the  total  mass  of  gas  in  the  upper  zone  (lbm) , 

11 
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r 
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Yn  is  the  ambient  oxygen  mass  fraction  (0.23), 

Ud 

mf  is  the  mass  flow  rate  of  fire  plume  gas  into 

the  upper  zone  (lbm/sec), 

is  the  total  oxygen  consumption  rate  by  all  fires 
(lbm/sec) , and 

mg  and  are  the  mass  flow  rates  of  gas  out  of  (ihg)  and 

into  (m^)  the  upper  zone  through  vents  (lbm/sec) . 

Following  the  development  given  in  [1]  for  the  upper  zone 
temperature,  the  oxygen  conservation  equation  can  be  simplified 
to  eliminate  ihg . Written  in  the  finite  difference  form  used 
in  the  computer  code,  the  expression  is 

(2-8)  AYq  = [(Y0a  - Yq)  (*f  + lhi)  - Q0]  t/MQ 

where  AYq  is  the  change  in  upper  zone  oxygen  concentration 
occurring  in  time  At. 

2.4  FORCED  VENTILATION 

Four  of  the  seven  full-scale  burn  tests  to  which  the  DACFIR 
model  was  compared  involved  forced  ventilation  of  the  cabin 
mock-up.  The  original  model  was  designed  to  account  for  gas 
flow  out  of  the  upper  zone  and  into  the  lower  zone  driven  only 
by  natural  buoyancy,  that  is  "natural"  ventilation.  Expressions 
were  adopted.  Equations  (5-12a)  and  5-12b)  of  [1] , using  the 
assumption  that  no  mechanically  imposed  pressure  difference 
exists  between  the  inside  and  outside  of  the  cabin.  An  attempt 
was  made  to  develop  expressions  that  would  account  for  both 
types  of  flows  simultaneously,  for  the  most  general  application, 
but  no  satisfactory  results  were  obtained.  DACFIR2  was  therefore 
revised  using  expressions  for  computing  the  flows  through  any 
vent  (doors,  escape  hatches,  or  ventilation  system)  based  totally 
on  the  pressure  difference  between  the  cabin  interior  and  exterior 
developed  by  the  ventilation  system  and  by  the  expansion  of  the 
upper  zone  gas.  A simple  orifice  flow  relationship  is  employed 
for  flow  through  a vent. 

(2-9)  *0  - 0.68Av(2APp/gc)1/2 
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where 


is  the  total  mass  outflow  rate  (lbm/sec) , 

2 

is  the  effective  area  of  the  vent  (ft  ) 


AP  = P . . - P . is  the  interior-exterior 

Cdbin  cunb  2 

pressure  difference  (lbf/ft  ) , 

g„  is  the  Newton  constant  (32.174  lbm-ft/lbf-sec2) , 

c 3 

and  p is  the  density  of  the  flowing  gas  (lbm/ft  ) . 

The  factor  0.68  is  an  orifice  coefficient.  For  a large  vent, 
such  as  a door,  the  effective  area  of  the  vent  is  determined  by 
computing  the  fraction  of  the  vent  within  either  the  upper  or 
lower  zone,  depending  upon  which  flow  is  being  calculated. 


Cabin  pressure  is  computed  from  the  gas  law  as 


(2-10) 


P , . = p RT 

cabin  u u 


where  pu  and  Tu  are  the  upper  zone  density  and  temperature  and 
the  gas  constant,  R,  has  the  value  for  air,  53.34  lbf-ft/lbm-°R. 

Expressions  of  the  form  of  Equation  (2-9)  are  used  for 
all  pressure  driven  flows  out  of  (or  into)  the  upper  and  lower 
zones,  the  densities  and  vent  areas  being  set  to  the  appropriate 
values  in  each  case.  In  the  case  of  specified  flow  rates  for 
certain  vents,  i.e.,  inflow  of  the  forced  ventilation,  the  mass 
flow  rate  produced  by  the  system  is  fixed  for  that  vent.  The 
ventilation  system  is  thus  assumed  to  be  able  to  overcome  the 
pressure  rise  caused  by  the  fire. 

2.5  THERMAL  DISCONTINUITY  POSITION 

The  original  DACFIR  model  was  designed  for  the  approximately 
rectangular  interior  of  wide-bodied  aircraft.  In  order  to 
simulate  fires  in  standard  width  cabin  structures,  the  computer 
code  was  modified  to  account  for  a more  cylindrical  geometry. 
Assuming  a truncated  circular  cabin  cross  section  (see  Figure  2-1) 
and  known  cabin  width  and  height,  expressions  were  derived  for 
the  cabin  volume  and  surface  area.  These  expressions  replaced 
their  rectangular-dependent  counterparts  in  the  equations  describing 
the  gas  dynamics.  In  particular,  a major  change  in  coding 


r 


resulted  from  the  thermal  discontinuity  position  (upper  layer 
depth)  being  implicitly  related  to  cross  sectional  area  rather 
than  by  the  explicit  relationship  of  the  original  model.  Accord- 
ingly, an  iterative  scheme  is  employed  to  solve  for  the  dis- 
continuity position. 

A switch  is  set  at  the  start  of  DACFIR2  indicating  selection 

of  rectangular  or  cylindrical  geometry.  The  equations  for  upper 

zone  volume,  V , and  upper  zone  surface  area,  A , remain 
u s u 

unchanged  from  the  original  model  in  the  rectangular  case.  In 
the  cylindrical  case,  the  upper  zone  volume  is 

(2-11)  Vu  = AxsCl  = [r2cos_1 (1-L/S) - (r-1) (2rL-L2) 1/2]  • CL 

where  L is  upper  zone  depth,  CL  is  cabin  length,  and  r is  a 
radius  defined  by 

r = (Cw2/4  + Ch2)/(2Cr) 

where  CH  is  the  cabin  height  and  Cw  the  width.  The  upper  surface 
area  is  given  by 

(2-12)  A = 2r  cos"1  (1-L/r)  • CT 

su  -u 

The  lower  zone  volume  is  calculated  as  the  total  volume 
less  the  upper  volume,  and  an  analogous  relationship  is  used 
to  obtain  the  lower  surface  area. 

Another  refinement  to  the  original  model  was  a correction 
for  the  seat  volume  contribution  to  the  total  cabin  volume 
which  is  a factor  in  the  upper  zone  depth  calculation.  The 
number  of  seat  groups  and  standard  seat  group  shape  are  used 
in  the  computation.  A series  of  tests  are  made  to  determine  the 
fraction  of  the  total  seat  volume  that  is  enveloped  by  the  upper 
zone  and  a correction  term  is  then  obtained  from  this  fraction 
for  both  the  upper  and  lower  zone  volumes.  With  the  corrected 
volumes  known,  the  thermal  discontinuity  is  obtained  from 
Equation  (2-11)  or  its  equivalent  for  the  rectangular  cross 
section  fuselage. 
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SECTION  3 


COMPARISON  OF  THE  MODEL 
TO  SEVEN  CABIN  FIRE  TESTS 


This  section  presents  the  results  of  the  comparison  of 

DACFIR2  to  seven  full-scale  cabin  mock-up  fire  tests  conducted 

during  the  Aerospace  Industries  Association  Crashworthiness 

r 4 51 

Development  Program  in  1968 1 ' . These  fire  tests  were  conducted 

to  determine  the  effect  of  furnishing  cabins  with  materials  of 
improved  fire  resistance.  Baseline  cases  were  established  by 
conducting  tests  with  cabins  furnished  with  materials  in  commercial 
use  at  the  time.  Several  sets  of  improved  materials  were  then 
tested  under  conditions  identical  to  the  baseline  tests.  A 
number  of  these  tests  involved  candidate  on-board  suppression 
systems  in  addition  to  improved  furnishing  materials . 

From  the  test  series,  six  tests  were  initially  selected  as 
most  appropriate  for  comparison  to  DACFIR2.  The  choice  was 
dictated  by  considerations  of  materials  used,  ventilation 
conditions,  and  whether  or  not  suppression  was  involved,  since 
the  present  DACFIR  model  does  not  include  fire  suppression.  A 
seventh  test  involving  only  an  ignition  source,  a pan  of 
jet  fuel,  was  added  later  as  a specific  exercise  for  the  gas 
dynamics  calculation.  After  the  validation  tests  (cases)  were 
selected,  the  furnishing  materials  involved  were  identified  and 
samples  of  the  materials  were  obtained  for  laboratory  testing. 

The  results  of  laboratory  tests  were  processed  into  a form 
compatible  for  input  to  the  DACFIR2  program.  Simulations  of 
each  case  were  made  and  the  results  put  into  tabular  and 
graphical  form  for  analysis. 

— 
f 4 1 

1 JNygren,  L.O.,  and  A.F.  Deardorff,  "AIA  Crashworthiness  Program  - 
Fire  Suppression  Section  - Results  of  Fire  Suppressant  System  Survey 
and  Tests",  D6-19456-3TN,  The  Boeing  Company,  November  1968. 

^Vaughn,  G.F.,  Deardorff,  A.F.,  and  L.O.  Nygren,  "AIA  Crashworthi- 
ness Program  - Fire  Suppression  Section  - Airplane  Crash  Fire 
Tests  to  Evaluate  Fire  Protection  Improvement  Devices  and  Materials", 
D6-19456-8TN,  The  Boeing  Company,  October  1968. 


3.1  THE  CABIN  MOCK-UP  FIRE  TESTS 


! 


The  seven  fire  tests  selected  for  validation  of  the  model 
may  be  divided  into  two  groups. 

3.1.1  15-Foot  Cabin  Mock-Up  Tests 

The  first  group  consists  of  five  tests  conducted  in 
a 15-foot  long,  constant  cross-section  stainless  steel  fuselage 
of  850  cubic  foot  volume.  Four  of  the  five  tests  involved 
furnishings.  The  fifth  test  - actually  the  averaged  results 
of  three  tests  - was  a fuel  pan  calibration  run  involving  only 
the  ignition  source  in  the  bare  metal  fuselage.  The  cabin 
interior  mock-up  in  the  15-foot  mock-up  tests  consisted  of  side- 
wall  panels,  a hatrack,  ceiling  panels  and  three  rows  of  simulated 
seats.  These  furnishings  were  installed  on  the  port  side  of  the 
fuselage  only,  the  starboard  side  and  floor  remained  bare  metal. 
Figure  3-1  shows  the  cross  section  of  the  15-foot  cabin  mock-up, 
and  Figure  3-2  gives  a plan  view  of  the  seat  arrangement. 

The  interior  furnishing  used  in  each  of  the  four 
15-foot  mock-up  tests  were  as  follows. 

Case  15P  - Present  In-Serivce  (1968)  Materials 

Sidewall:  Vinyl-coated  aluminum  sheet. 

Hatrack:  Polyurethane  foam  covered  with 

vinyl-coated  fiberglass  cloth. 

Ceiling:  Paper-honeycomb-core,  fire-retardant 

polyester-f iberglass-laminate-faced 
sandwich  panel  covered  with  vinyl- 
coated  fiberglass. 

Seats:  Muslin-covered  polyurethane  foam, 

upholstered  in  Nylon  fabric. 


Case  15A  - Improved  Materials  Set  "A" 

Sidewall:  Enameled  aluminum  panels  with 

PVF*  surface  film.  (Samples 
of  this  material  were  not  avail- 
able for  laboratory  testing) . 

•Polyvinyl  Fluoride 
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Hatrack : 


Ceiling: 


Seats : 


Nylon-honeycomb-core , fire- 
retardant  epoxy  f iberglass-laminate- 
faced  sandwich  panel  covered  with 
a fiberglass  mat  for  padding  and  a 
final  covering  of  vinyl-coated 
fiberglass . 

Nylon-honeycomb-core , fire-retardant 
epoxy  fiberglass  laminate-faced 
sandwich  panel  with  printed  PVF 
surface. 

Muslin-covered  polyurethane  foam 
upholstered  with  Nylon  fabric. 


Case  15B  - Improved  Materials  Set  "B" 


Sidewall  and  Ceiling: 


Hatrack: 


Seats : 


Nylon-honeycomb-core , fire-retardant 
epoxy  fiberglass-laminate-faced 
sandwich  panel  with  printed  PVF 
surface. 

Nylon-honeycomb-core , fire-retardant 
epoxy  fiberglass-laminate-faced 
sandwich  panel  with  PVF. 

Flame-retardant  polyurethane 
covered  with  a Nylon-fabric  slip- 
cover, and  upholstered  with  Nylon 
fabric. 


Case  15C  - Improved  Materials  Set  "C" 

Sidewall  and  Ceiling:  Nylon-honeycomb-core,  fire-retardant 

epoxy  fiberglass-laminate-faced 
sandwich  panel  with  printed  PVF 
surface. 

Hatrack:  Nylon-honeycomb-core,  fire-retardant 

epoxy  fiberglass-laminate-faced 
sandwich  panel  with  PVF  surface. 

Seats:  Fiberglass  cushion  in  a Nylon  slip- 

cover, upholstered  with  Nylon 
fabric. 

Materials  in  the  first  case,  15P,  were  characteristic 
of  the  then  in-service  materials  which  involved  extensive  use  of 
vinyl  coated  fabrics,  paper  honeycomb  structural  panels,  and 
non-fire-retardant  polyurethane  (PU)  foams.  The  improved  materials 
of  Case  15A  included  hatracks  and  ceilings  of  a nylon  based  honey- 
comb sandwich.  The  sidewall  in  this  case  was  of  an  enameled 
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aluminum  sheet  construction,  samples  of  which  were  not  available 
for  laboratory  testing.  Lack  of  this  data  may  not  completely 
invalidate  comparison  of  the  program  to  this  case  for  reasons 
mentioned  below.  The  materials  of  Cases  15B  and  15C  were  identical 
except  for  the  seat  cushion  material:  a fire  retardant  poly- 
urethane foam  in  15B  and  an  inert  fiberglass  padding  in  15C. 
Comparison  of  the  results  of  these  two  tests  indicate  the  effect 
of  the  seat  foam  material  which  was  suspected  of  being  a major 
contributor  to  the  fire. 

Ventilation  of  the  15-foot  cabin  mock-up  was  provided 
by  the  forced  flow  from  a regulated  blower.  Air  entered  the 
mock-up  at  a baffeled  inlet  duct  forward  and  exited  through  a 
vent  near  ceiling  level  at  the  rear.  A flow  rate  of  200  cubic 
feet  per  minute  was  maintained  by  the  blower  throughout  the 
tests . 

Ignition  of  the  fire  in  each  of  the  five  15-foot 

mock-up  tests  was  accomplished  in  the  same  manner,  a 12  inch  by 
2 

12  inch  (1  ft  ) metal  pan  filled  with  one  quart  of  Jet  A-l 
fuel  floating  on  water  was  placed  beneath  the  center  seat  row 
next  to  the  sidewall  (see  Figure  3-2) . From  10  to  30  seconds 
were  required  for  flames  to  spread  over  the  entire  fuel  pool 
surface  after  ignition  was  initiated  at  one  edge  of  the  pan. 

When  the  fuel  surface  was  completely  involved  in  flame,  that  time 
was  denoted  as  zero  in  the  data  collection  process.  The  fire 
in  each  test  was  allowed  to  burn  for  300  seconds  (five  minutes) 
after  time  zero. 

Instrumentation  in  each  of  the  15-foot  mock-up  tests 
provided  measurements  of  cabin  atmosphere  temperature  and 
composition.  Temperature  was  measured  by  six  thermocouples: 
three  located  two  inches  below  the  ceiling  on  the  cabin  fore- 
aft  centerline  at  one  inch,  99  inches,  and  144  inches  from  the 
aft  bulkhead;  and  three  located  two  inches  below  the  left 
hatrack  at  43  inches,  99  inches,  and  133  inches  from  the  aft 
bulkhead.  Figure  3-2  shows  the  position  of  the  six  thermocouples. 
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The  gas  sampling  port  was  located  at  42  inches  above 
the  floor  and  the  seine  distance  forward  of  the  aft  bulkhead.  The 
end  of  the  tube  forming  the  port  extended  approximately  24 
inches  into  the  cabin  from  the  starboard  sidewall  (see  Figures 
3-1  and  3-2) . Gas  samples  were  withdrawn  through  the  port  and 
collected  in  plastic  bags  for  chemical  analysis.  The  cabin 
atmosphere  was  sampled  for  oxygen,  carbon  dioxide,  carbon  monoxide, 
hydrocarbons,  acid  gases,  chlorides,  fluorides,  acrolein. 

Freon  1301,  ammonia,  nitrogen  dioxide,  carbonyls,  and  hydrogen 
cyanide.  Concentrations  of  these  gases  were  reported  at  30 
second  intervals.  The  laboratory  tests  on  the  furnishing  materials 
measured  only  the  combustion  products  carbon  dioxide,  carbon 
monoxide,  hydrogen  chloride,  hydrogen  fluoride,  hydrogen 
cyanide,  and  nitrogen  dioxide  and  thus  only  the  concentrations 
of  these  species,  in  addition  to  oxygen  concentration,  were  used 
in  comparisons  to  the  computer  simulations.  Table  3-1  gives  the 
measurement  technique  in  the  full  scale  testing  for  each  sub- 
stance and  the  estimate  of  measurement  accuracy. 

TABLE  3-1 

ANALYSIS  OF  THE  CABIN  ATMOSPHERE  - AIA  CABIN  FIRE  TESTS 


Substance 


Method  of  Analysis 


Accuracy 


Oxygen 

Carbon  Dioxide 
Carbon  Monoxide 
Hydrogen  Chloride 
Hydrogen  Fluoride 

Nitrogen  Dioxide 
Hydrogen  Cyanide 


Potentiometric  Oxygen  Meter 

Volumetric  Meter 

Detector  tubes 

Silver  nitrate  titration 

Thorium  nitrate-alzarin 
titration 

Detector  tubes 

Detector  tubes 


+ 0.5%  02 
+ 0.2%  by  volume 

+ 25%  by  volume 
* 

* 

+ 25%  by  volume 
+ 25%  by  volume 


* 

accuracy  estimates  for  these  substances  not  given  in  terms 
of  a percentage  of  the  measurement. 
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3.1.2  26-Foot  Cabin  Section  Tests 


The  remaining  two  tests  of  the  seven  selected  for  com- 
parison to  the  model  were  conducted  in  the  aft  26  feet  of  a 
B727  fuselage  cabin  furnished  to  simulate  a complete  passenger 
cabin.  Total  cabin  volume  was  about  2000  cubic  feet.  The  first 
of  these  two  tests,  designated  26P,  used  furnishing  characteristic 
of  the  then  in-service  materials.  The  other  test,  26N,  employed 
a set  of  improved  materials.  Materials  used  in  the  26  foot 
fuselage  cases  were  similar  to  those  in  Cases  15P  and  15B  with 
some  additions  such  as  carpet  and  passenger  service  units. 

The  materials  for  each  test  were  as  follows. 


Case  26P  - Present  In-Service  (1968)  Materials 


Carpet: 


Sidewall : 

Window  Reveals : 
Window  Dust  Panes: 
Hatrack: 


Passenger  Service  Units: 
Ceiling  Panels: 

Seats: 


Acrilan  pile  with  muslin  faced 
polyurethane  foam  pad. 

Vinyl-Aluminum  laminate. 


Acrylic 

Vinyl  coated  fiberglass,  poly- 
urethane foam,  resin  coated 
fiberglass,  and  resin  impregnated 
paper  honeycomb.  Bull  nose:  ABS 
covering  polyurethane  foam. 


Paper  honeycomb  core,  resin 
coated  fiberglass  fabric,  vinyl 
cover . 

Wool  and  muslin  upholstered 
polyurethane  foam. 


Case  26N  - Improved  Materials 


Carpet: 


Sidewall: 


Window  Reveals: 


High  temperature  resistant 
nylon  pile  (material  not  available 
for  laboratory  testing) . 

PVF  faced  nylon  honeycomb  core 
sandwich  panels. 

Polycarbonate 


Window  Dust  Panes: 


Polycarbonate  and  acrylic 
(alternating) 

Hatrack:  PVF  faced  nylon  honeycomb  core 

sandwich  panels. 

Passenger  Service  Units:  Polycarbonate  (units  simulated 

by  sheet  stock) . 

Ceiling  Panels:  PVF  faced  nylon  honeycomb  core 

sandwich  panels. 

Seats:  High  temperature  resistant  nylon 

and  muslin  upholstered  poly- 
urethane foam. 

Figure  3-3  shows  a cross-section  of  the  727  cabin  and  Figure  3-4 
gives  a plan  view  of  the  seat  locations  for  both  tests. 

Ventilation  of  the  cabin  was  by  natural  flow  out 
two  openings,  a 2 foot  wide  by  4 foot  high  exit  at  the  aft 
and  the  second  a 3 foot  wide  by  7 foot  opening  in  the  forward 
bulkhead.  The  ignition  source  was  a 30  x 30  inch  fuel  pan 
containing  seven  quarts  of  Jet  A-l  fuel  on  water.  An  18  by 
24  inch  section  of  the  cabin  wall  was  cut  out  at  floor  level 
and  a section  18  inches  wide  and  8 inches  deep  removed  from 
the  floor  to  simulate  a fuselage  rupture.  The  fuel  pan  was 
inserted  in  this  region  extending  about  half  way  into  the 
cabin.  Figure  3-5  shows  this  arrangement.  Time  zero  for 
measurements  in  the  26-foot  cabin  tests  was  intended  to  be,  as 
in  the  15-foot  mock-up  cases,  the  time  at  which  the  entire  fuel 
pan  was  involved  in  flames.  This  procedure  was  followed  for 
26P  but  problems  were  encountered  with  Case  26N.  In  this  case 
flames  would  not  enter  the  cabin  but  remained  totally  on  the 
outside  of  the  fuselage  until  an  exterior  wind  screen  was 
removed  at  78  seconds  after  ignition.  At  that  time  flames  did 
enter  the  cabin  and  came  in  contact  with  the  interior  materials. 
Time  calibration  of  data  in  case  26N  was  therefore  adjusted  so 
that  time  zero  was  the  point  when  flames  entered  the  cabin.  To 
prevent  damage  to  the  fuselage  structure,  case  26P  was 
terminated  with  CC^  extinguishment  at  202  seconds  and  case  26N 
at  162  seconds  (adjusted) . 
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Figure  3-3.  Cross  Section  of  the  26  Foot  Fuselage  Test  Section 
Only  the  position  of  the  port  seat  row  is  shown. 


•Area  modeled  by  DACFIR2  for  fire  spread 


Figure  3-4.  Plan  View  of  the  26  Foot  Fuselage  Test  Section 


More  instrumentation  was  used  in  the  26-foot  fuselage 
tests  than  in  the  15- foot  mock-up  tests.  The  temperature  field 
of  the  cabin  atmosphere  was  measured  by  50  thermocouples  located 
on  five  trees.  Each  tree  provided  measurements  at  the  10  points 
shown  in  Figure  3-6.  The  position  of  each  tree  in  the  cabin  is 
shown  in  Figure  3-4.  Gas  sampling  was  made  through  two  ports  in 
the  furnished  section.  Port  A was  48  inches  from  the  floor  and 
directly  opposite  the  ignition  fire  and  port  B was  at  the  same 
vertical  level  and  on  the  same  side  as  A but  near  the  forward 
bulkhead.  These  positions  are  shown  in  Figure  3-4.  Samples 
of  the  cabin  atmosphere  were  withdrawn  at  30  second  intervals 
and  analyzed  for  the  same  compounds  as  the  samples  in  the  15- 
foot  mock-up  cases  (see  Table  3-1) . Additional  temperature 
measurements  were  made  in  the  26-foot  fuselage  tests  using  shielded 
and  aspirated  thermocouples,  calorimeters,  and  wet  bulb 
thermocouples . 

3.2  LABORATORY  TESTS  OF  CABIN  MATERIALS 

Furnishing  materials  of  the  AIA-CDP  cabin  mock-up  tests 
described  above  were  subjected  to  a series  of  laboratory  tests 
to  provide  input  data  for  the  DACFIR  computer  program.  Tests 
were  conducted  under  subcontract  to  UDRI  by  the  Boeing  Materials 
Services  division  of  the  Boeing  Commerical  Airplane  Company. 

The  final  technical  report  submitted  by  Boeing  is  attached  as 
Appendix  B of  this  report  and  contains  a full  description  of 
the  testing  program.  Testing  was  conducted  on  samples  of  25 
separate  materials.  Three  other  materials  used  in  the  full 
scale  burn  tests  could  not  be  obtained  for  laboratory  analysis. 

Tests  were  made  in  the  Ohio  State  University  (OSU)  combustion 
analyzer  for  the  following: 

o Rate  of  surface  flame  spread  in  horizontal,  vertically 
upward,  and  vertically  downward  directions; 

o Time  to  flame  (ignite)  when  exposed  to  a small 
pilot  flame; 
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Figure  3-5.  Fuel  Pan  Location  Details  - 
26-Foot  Cabin  Cases 


Figure  3-6. 


Thermocouple  Locations  in  the  26-Foot 
Fuselage  Test  Section 
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o Time  to  burn  (time  from  ignition  to  the  burned  out  or 
charred  state) ; 

o Heat  release  rate  and  total  heat  released; 
o Smoke  release  rate  and  total  smoke  released. 

Each  of  these  items  was  measured  at  four  levels  of  externally 

2 2 2 

applied  heat  flux,  1.32  Btu/ft  -sec  (1.5  W/cm  ),  2.20  Btu/ft  -sec 

(2.5  W/cm2).  3.08  Btu/ft2-sec  (3.5  W/cm2),  and  3.96  or  4.40 
0 2 

Btu/ft  -sec  (4.5  or  5.0  W/cm  ).  Materials  were  tested  in  the 
OSU  analyzer  in  the  orientation (s)  in  which  they  were  installed 
in  the  cabin. 

Tests  were  also  conducted  on  the  materials  in  the  smoldering 
state  (radiant  exposure  only  without  pilot  flame)  in  the  OSU 
device  for  the  data  items. 

o Time  to  smolder  (time  from  start  of  exposure  at  given 
flux  level  to  the  start  of  pyrolysis) 

o Time  to  char  when  smoldering  (time  from  the  start  of 
smoldering  to  the  finish) 

o Rate  of  smoke  release  in  the  smoldering  state. 

Rate  of  heat  release  in  the  smoldering  state  was  also  measured 
but  this  quantity  is  not  used  by  the  present  DACFIR  model. 
(Examination  of  the  measured  rates  of  smoldering  heat  release 
show  that  they  are  almost  always  insignificant  with  respect  to 
the  flaming  heat  release  rates,  thus  justifying  their  neglect.) 

Flame  spread  rates  and  the  various  times  of  transition 
were  recorded  by  Boeing  and  reported  for  individual  runs  and 
averages  of  several  runs  when  more  than  one  run  was  made.  Rate 
of  heat  release  and  rate  of  smoke  release  data  from  the  OSU 
analyzer  were  recorded  on  strip  charts.  The  curves  from  these 
charts  were  automatically  digitized  and  input  to  a data  processing 
program  that  determined  the  total  amount  of  heat  and  smoke 
release  (by  integration)  and  the  maxima  and  time  averages  of  these 
quantities.  Results  from  multiple  runs  were  also  averaged  and 
reported. 
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Testing  of  the  interior  materials  for  toxix  gas  generation 
was  done  using  a National  Bureau  of  Standards  (NBS)  Smoke 
Chamber  specially  modified  for  gas  composition  measurements. 
Concentrations  of  HC1,  HF,  HCN,  N02,  and  S02  were  measured 
with  colorimetric  indicator  tubes  while  C02  and  CO  were 
measured  by  process  instruments.  Appendix  B gives  the  details 
of  these  measurements. 

Additional  processing  of  the  materials  data  supplied  by 
Boeing  was  necessary  in  order  to  form  input  data  for  the  DACFIR 
model.  The  type  of  treatment  differed  for  each  data  item  and 
is  best  described  by  considering  each  item  in  turn. 

Flame  spread  rates  Flame  spread  rates  were  converted  to  units 
of  feet/sec  from  the  reported  units  of  inches/min. 


Time  to  flame.  Time  to  burn  Two  methods  were  used  to 
determine  the  times  of  ignition  and  burnout.  First  a visual 
observation  of  these  times  was  made  for  each  test.  Second, 
the  heat  release  rate  curve  was  examined  to  determine  when, 
for  example,  a net  positive  rate  of  heat  release  indicated 
ignition  or  when  the  rate  of  heat  release  returned  to  zero 
indicating  burn-out.  While  the  second  method  appears  at  first 
to  be  more  accurate,  it  was  found  that  the  thermal  inertia  of 
the  OSU  analyzer  was  so  significant  that  times  read  from  the 
raw  heat  release  curves  could  lag  behind  the  visual  measurements 
as  much  as  30  seconds  or  more.  The  visual  measurements  were 


therefore  used  in  all  cases. 

Heat  Release  Rate  The  Boeing  computer  program  for  processing 
the  OSU  analyzer  data  integrates  the  heat  release  rate  curve 

to  obtain  total  heat  release.  This  value  is  reported  in 

2 2 
Joules/cm  . After  conversion  to  units  to  Btu/ft  , the  heat 

release  rate  for  each  material  was  computed  by  dividing  the 

total  heat  released  by  the  time  to  burn  (average  of  all  values 

if  more  than  one  run  was  made) . The  heat  release  rate  used 

is  therefore  an  average  value,  constant  over  the  sample  burn 

time. 
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Smoke  Release  Rate  Cumulative  smoke  release  was  reported  as 
optical  density  per  unit  sample  area.  This  number  was  con- 
verted to  "particles"  of  smoke  released  per  unit  area  by 

2 

dividing  by  the  conversion  factor  0.045757  ft  /particle.  An 
average  release  rate  was  then  found  by  dividing  the  total 
smoke  particles  released  per  unit  area  by  the  time  to  burn 
or  time  to  smolder  out  as  appropriate. 

Toxic  Gas  Release  Rates  The  NBS  Smoke  Chamber  is  not  a flow- 
through device;  instead  smoke  and  gases  emitted  by  the  sample 
accumulate  within  the  18-cubic  feet  of  the  chamber  volume. 

Since  a gas  release  rate  is  to  be  supplied  as  input  to  DACFIR, 
several  operations  were  required  to  convert  the  reported  con- 
centration values  to  release  rates.  Concentrations  of  each  gas 
specie  of  interest  were  measured  at  1,  2,  4,  and  10  minutes  during 
the  tests.  The  maximum  concentration,  usually  that  at  10  minutes, 
was  converted  to  mass  units  using  the  chamber  volume  and 
assuming  the  chamber  gas  to  be  at  atmospheric  pressure  and 
room  temperature.  The  rate  of  emission  per  unit  sample  area 
was  then  computed  from  the  reported  value  of  time  to  burn  out 
or  time  to  smolder  out  and  the  sample  area.  Corrections  to 
the  release  rates  were  made  for  two  gases,  CC>2  and  CO. 

When  samples  were  tested  in  the  flaming  mode,  C02  was 
generated  by  the  propane  pilot  flame.  This  rate,  known 
independently  from  calibration  runs  using  an  inert  sample, 
is  about  340  ppm  per  minute.  A small  correction  was  also 
made  for  the  atmospheric  concentration  of  C02 , 300  ppm.  Gen- 
eration of  carbon  monoxide  by  the  pilot  flame  was  also 
detected  in  calibration  runs.  The  rate  of  generation  varied 
with  the  radiant  exposure  level,  decreasing  at  higher  fluxes. 
(Indicating,  perhaps,  that  CO  is  produced  by  quenching  of  the 

pilot  on  the  cool  sample  surface) . The  maximum  rate  of  150 

2 2 

ppm/min  was  detected  at  a flux  of  1.32  Btu/ft  -sec  (1.5  W/cm  ) 

2 2 

and  the  smallest  rate  of  30  ppm  at  4.4  Btu/ft  -sec  (5.0  W/cm  ). 
These  rates  were  used  to  correct  all  CO  emission  data  for 
flaming  exposures. 
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3.3  COMPARISON  OF  THE  DACFIR  MODEL  TO  THE  TEST  RESULTS 

The  results  of  the  DACFIR2  model's  simulation  of  each  of  the 
seven  cabin  fire  tests  are  discussed  below.  The  comparison  of 
the  program  output  to  the  test  data  is  given  in  two  forms : 

(1)  plots  of  the  predicted  and  measured  values  of  cabin  gas 
temperature  and  composition;  and  (2)  tables  and  figures  giving 
the  areas  of  burned  material  computed  by  the  model  and  observed 
at  the  end  of  each  test.  Figures  are  also  presented  for  smoke 
concentration,  but  since  no  quantitative  smoke  measurements  were 
made  in  the  tests  only  the  model's  results  are  shown.  The 
results  presented  are  those  after  all  the  refinements  described 
in  Section  2 were  made  and  represent  the  current  ability  of 
DACFIR2  to  reproduce  these  tests  based  upon  the  material  input 
data  used. 

3.3.1  Case  15Z  - Fuel  Pan  Calibration  Tests 

The  first  comparison  of  model  and  tests  was  for 
a series  of  "fuel  pan  calibration"  tests  conducted  in  the  15- 
foot  mock-up  fuselage  without  furnishing  materials.  The  fire 
consisted  of  the  one  square-foot  pan  of  Jet-A  fuel.  This  fuel, 
fire  size,  and  location  were  subsequently  used  as  the  ignition 
source  in  each  test  with  materials.  The  ability  of  the  DACFIR 
program  to  correctly  simulate  these  "no-spread"  burn  tests  was 
used  as  an  intermediate  step  in  validating  the  program. 

Heat  and  smoke  release  rates  for  jet  fuel  were 
measured  in  the  OSU  analyzer  and  the  values  were  used  as  input 
for  the  simulation  of  the  fuel  pan  fire.  Data  on  Jet  A from 
Sarkos ^ was  used  to  estimate  the  carbon  monoxide  production 
rate  for  the  fuel.  Since  no  direct  measurements  of  C02  produc- 
tion could  be  found  for  jet  fuel,  this  rate  was  estimated  by 
assuming  a simple  combustion  model  (reaction  equation)  for  the 
jet  fuel.  The  rate  of  release  of  carbon  as  smoke  was  calculated 


Sarkos,  C.P.  "Measurement  of  Toxic  Gases  and  Smoke  from 
Aircraft  Cabin  Interior  Materials  Using  the  NBS  Smoke  Chamber 
and  Colorimetric  Tubes",  FAA-RD-76-7,  March  1976. 
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from  the  observed  smoke  release  rate  by  a method  of  Seader  and 
f 8 1 

0ul  . With  the  rate  of  carbon  release  known  for  CO  and  smoke, 
the  remaining  carbon,  about  94%  of  the  total  in  the  fuel,  was 
assumed  to  occur  as  CC^-  Since  the  mass  burning  rate  for  the 
fuel  is  specified  as  input  for  the  model  a C02  rate  could 
then  be  computed. 

Figure  3-7  shows  the  computed  upper  zone  temperature, 
solid  curve,  and  the  readings  from  two  representative  ceiling 
thermocouples,  the  broken  curves.  The  pattern  of  temperature  rise 
in  the  test  is  a general  increase  from  the  ambient,  about  70 °F, 
to  between  250°F  at  the  aft  ceiling  and  340°F  at  the  center 
ceiling.  The  theoretical  curve  rises  quickly  to  218 °F  and  then 
maintains  a very  slow  rise  to  245 °F  when  the  test  was  terminated. 
Figures  3-8,  3-9,  and  3-10  show,  respectively,  the  comparisons  of 
measured  and  computed  oxygen,  carbon  dioxide,  and  carbon 
monoxide  concentrations.  In  each  figure  the  observed  and 
computed  results  seem  to  follow  the  same  trends  of  slow  increase, 
or  decrease  for  0 2,  throughout  the  run.  The  computed  upper 
zone  depth  of  this  run  (not  shown)  was  quickly  established 
at  about  3.5  feet  and  grew  slowly  to  4.5  feet  by  300  seconds. 

The  figures  show  generally  good  agreement  between 
the  model's  predictions  and  the  measured  results  for  this 
fairly  simple  case  in  which  there  is  no  fire  spread.  In  the 
comparison  of  temperatures.  Figure  3-7,  the  rapid  rise  of  the 
computed  temperature  is  due  to  the  assumptions  of  a steady- 
state  fire  and  plume  and  a spatially  uniform  upper  zone  which 
does  not  mix  with  the  lower  zone.  Most  of  the  difference 
between  the  computed  and  observed  results  for  the  first  two 
minutes  is  probably  due  to  these  assumptions  since  in  the 
test  the  early  gas  motion  in  the  upper  cabin  would  be  expected 
to  include  significant  mixing.  The  approximately  constant 
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1 J Seader,  J.D.,  and  S.S.  Ou,  "Correlation  of  the  Smoking 
Tendency  of  Materials,"  Fire  Research,  Vol.  1,  No.  1,  March  1977. 


difference,  about  75°F,  between  the  aft  and  center  ceiling 
temperatures  throughout  the  test  indicates  that  some  cooling 
of  the  horizontal  ceiling  flow  occurs  even  in  this  short  fuselage. 

The  model  did  well  in  predicting  the  accumulation 
of  carbon  monoxide  and  the  consumption  of  oxygen  and  somewhat 
less  well  the  C02  concentration.  Since  the  method  of  establishing 
the  input  value  of  C02  production  rate  was  rather  approximate, 
the  difference  in  the  results  is  not  surprising. 

3.3.2  Case  15P  - Present  In-Service  (1968)  Materials 

Figures  3-11  through  3-20  give  the  calculated  and 
observed  results  from  Case  15P,  the  pre-1968  furnishings  in 
the  15  foot  mock-up  cabin.  Figure  3-11  shows  that  the  model 
predicts  the  upper  cabin  temperature  reasonably  well  for  the 
first  three  minutes  of  the  test.  After  180  seconds  there  seems 
to  have  been  a flashover  in  the  cabin  mock-up  as  indicated  by 
the  sharp  rise  of  the  themocouple  readings  at  200  seconds. 

After  the  flashover  the  fire  diminished  substantially  or  self- 
extinguished  as  indicated  by  the  falling  temperature  readings. 

The  model  does  not  consider  flashover  (the  mechanism  is  as  yet 
not  well  understood)  and  so  the  computed  temperature  rises 
slowly  toward  an  equilibrium  of  about  1050 °F  during  the  last 
two  minutes  of  the  test.  It  is  interesting  to  note,  however, 

that  if  one  simple  criterion  for  flashover  is  used,  that  given 

r 2 1 

by  Quintiere1  as  an  upper  zone  temperature  of  600°C  (1112°F) , 
the  gas  temperature  predicted  by  DACFIR2  reached  the  vicinity 
of  1050 °F  at  about  215  seconds  and  thus  the  model  indicates 
that  flashover  could  have  been  expected  at  or  near  200 
seconds . 

In  Figure  3-12  the  computed  and  observed  oxygen 
concentrations  are  compared.  It  can  be  seen  that  the  computed 
curve  predicts  a much  lower  oxygen  concentration  than  that 
ovserved.  The  discrepancy  is  probably  due  to  the  assumption 
of  complete  (stoichiometric)  burning  by  the  model  and  so  suggests 
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that  the  actual  combustion  efficiency  is  rather  low.  For  the 
toxic  gas  accumulation  Figures  3-13  through  3-17  show  that, 
with  the  exception  of  HC1,  the  concentrations  are  reasonably 
well  predicted  throughout  the  test.  Once  again  carbon  dioxide 
is  overestimated,  and  this  could  be  primarily  the  effect  of  the 
ignition  jet  fuel  fire.  Figure  3-18  shows  the  predicted  upper 
(hot  gas)  zone  depth  and  Figure  3-19  the  predicted  smoke  accumu- 
lation. As  previously  mentioned,  no  test  results  were  available 
for  comparison  for  these  items.  The  figures  show  that  the  entire 
cabin  volume  fills  with  hot  gas  at  about  110  seconds  and  that 
at  about  the  same  time  the  light  transmission  over  one  foot 
drops  below  one  percent. 

Figure  3-20  shows  the  predicted  spread  of  the  areas 
of  burning  material  at  several  times  in  the  simulation.  The 
shaded  region  in  the  same  figure  shows  the  observed  fire  damage 
at  the  end  of  the  test.  Table  3-2  gives  a quantitative  compar- 
ison of  the  involvement  of  each  surface  as  determined  from  the 
post-test  damage  assessment  and  by  the  progress  of  the  fire  in 
the  simulation.  Only  the  section  used  in  the  model  is  considered 
in  the  figure  and  table.  This  section  normally  contains  all 
the  floor,  seat,  and  sidewall  damage  in  the  full  scale  test, 
while  ha track  and  ceiling  damage  can  extend  out  of  the  section. 
The  figure  and  table  show  that  the  model  predicts  well  the  total 
involvement  of  the  hatrack  and  ceiling  but  generally  overesti- 
mates the  spread  on  the  sidewall  and  seats.  No  spread  could, 
of  course,  occur  on  the  inert  floor.  The  overestimate  of  flame 
spread  on  the  chlorine  containing  vinyl  sidewall  is  probably 
the  reason  (barring  measurement  error)  for  the  disagreement  in 
the  HC1  prediction.  The  figure  also  shows  the  flame  spread 
from  seat  to  seat  in  the  row  exposed  to  the  ignition  source. 

The  model  predicts  an  accelerating  spread  of  flames  over  the 
seat  row,  while  in  the  test  the  seat  fire  was  more  confined, 
involving  mainly  the  seat  over  the  fuel  pan  and  the  adjacent 
seat. 
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Figure  3-13.  Carbon  Dioxide  Concentration  - Case  15P 
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Figure  3-16.  Hydrogen  Cyanide  Concentration  - Case  15P 
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Figure  3-19.  Smoke  Accumulation  - Case  15P 


TABLE  3-2 


COMPARISON  OF 
OF 

SIMULATED  AND 
DAMAGE  - CASE 

OBSERVED  AREAS 
15P 

OBSERVED 

AREA  ^ 1 ^ 

COMPUTED 

AREA 

STRUCTURE 

(ft2) 

% TOTAL 
EXPOSED 

(ft2) 

% TOTAL 
EXPOSED 

Sidewall 

17  to  26 

20  to  30 

45.0 

100 

Ha track 

32.5 

100 

33.8 

100 

Ceiling 

59.0 

100 

86.0(2] 

100 

Seats 

25.0 

55 

49.5 

100 

Damaged  area  within  test  section  considered  by  model 
(center  7.5  ft  of  15  ft  length) 

r ? I 

Difference  in  percent  figures  due  to  differences  in 
modeled  and  actual  ceiling  area 


3.3.3  Case  15A  - Improved  Materials  Set  A 

Figure  3-21  shows  the  predicted  and  observed  cabin 
temperature  for  this  test  using  improved  ceiling  and  hatrack 
materials.  Since  the  experimental  PVF-aluminum  material  was 
not  available  for  laboratory  testing,  an  inert  material 
was  specified  for  the  sidewall  in  the  computer  simulation. 

Without  the  active  involvement  of  the  sidewall,  the  fire  as 
computed  by  DACFIR  was  confined  to  the  center  seat  group,  the 
hatrack,  and  the  ceiling.  Figures  3-22  through  3-27  show  the 
comparison  of  the  oxygen  and  toxic  gas  levels.  The  agreement 
for  02,  CC >2,  and  HCN  is  good  for  most  times,  while  that  for  CO 
and  HC1  is  fair  or  poor  before  three  minutes  but  improves  after 
that  time.  HF  is  not  predicted  well  during  the  test  but  this  is 
to  be  expected  since  the  PVF  coated  sidewall  data  was  not  avail- 
able. The  computed  smoke  accumulation  and  the  growth  of  the  hot 
layer  differed  little  from  the  results  for  Case  15P  as  can  be 
seen  by  comparing  Figures  3-28  and  3-29  to  Figures  3-18  and  3-19. 

Figure  3-30  shows  the  observed  and  computed  damage 
sustained  in  Case  15A.  The  computed  results  are  confined  to 
the  lining  surfaces  of  the  hatrack  and  ceiling  by  the  inert 
floor  and  sidewall.  Table  3-3  compares  the  areas  of  damage  in 
absolute  terms  and  as  percentages  of  each  surface  exposed. 

Without  the  fire  spread  to  the  forward  and  aft  seat  groups  from 
the  sidewall,  a closer  agreement  in  the  total  involvement  of 
seat  material  was  obtained  than  in  Case  15P. 

3.3.4  Case  15B  - Improved  Materials  Set  B 

In  this  test  the  seat  upholstery  and  cushioning  were 
of  newer,  more  fire-resistant  materials.  In  addition,  the 
sidewall  was  made  of  the  same  nylon-honeycomb  sandwich  panel 
construction  as  the  ceiling  and  hatrack  in  this  and  the  previous 
improved  materials  test. 

Figure  3-31  shows  the  computed  and  observed 
temperature.  Here  the  model  did  not  reproduce  the  general  rise 
of  temperature  starting  at  about  100  seconds.  The  observed 


41 


/ 


45 


r 


TABLE  3-3 

COMPARISON  OF  SIMULATED  AND  OBSERVED  AREAS 
OF  DAMAGE  - CASE  15A 

OBSERVED  AREA  ^ 1 ^ COMPUTED  AREA 


STRUCTURE 

Hi 

ft 

N> 

% TOTAL 
EXPOSED 

(ft2) 

% TOTAL 
EXPOSED 

Sidewall 

34.6 

40 

Material  Data  Not 
Available 

Ha track 

32.5 

100 

33.8 

100 

Ceiling 

59.0 

100 

86.0^ 

99 

Seats 

22.4 

so'11 

44.5 

90 

11 J Damaged  area  within  test  section  considered  by  model 
(center  7.5  ft  of  15  ft  length) 

^ Difference  in  percent  figures  due  to  differences  in 
modeled  and  actual  ceiling  area 
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Figure  3-30.  Computed  and  Observed  Areas  of  Damage  - Case  15A 
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Figure  3-31.  Upper  Zone  Temperature  - Case  15B 


temperature  stabilizes  after  about  200  seconds  and  only  near 
the  end  of  the  test  does  the  computed  temperature  approach  the 
measured  levels.  The  oxygen  and  CC>2  concentrations  shown  in 
Figures  3-32  and  3-33  are  also  not  in  as  good  agreement  as  in 
Cases  15P  and  15A.  For  carbon  monoxide.  Figure  3-34,  the 
agreement  is  better,  equaling  that  obtained  in  15P  and  15A. 

The  match  for  HC1  is  difficult  to  judge.  No  HC1  was  detected 
before  210  seconds  but  after  this  time  the  measured  values  jump 
to  over  14,000  ppm.  Somewhat  the  same  rapid  variation  in  the 
measured  values  of  HCN  and  particularly  HF  is  shown  in  Figures 
3-36  and  3-37.  The  agreement  with  the  computed  values  can  only 
be  judged,  at  best,  fair.  The  time  required  for  the  cabin  to 
fill  with  the  combustion  products  was  somewhat  longer,  135 
seconds  versus  about  100  seconds  in  the  earlier  tests,  as 
shown  in  Figure  3-38.  Smoke  accumulation  was  also  slower,  the 
visibility  dropping  below  one  percent  at  150  seconds. 

In  Table  3-4  and  Figure  3-40  the  pattern  of  good 
prediction  of  involvement  on  the  ceiling  and  hatrack  and  over 
prediction  on  the  sidewall  and  seats  is  again  evident. 

3.3.5  Case  15C  - Improved  Materials  Set  C 

The  third  set  of  improved  materials  consisted  of  those 
used  in  Case  15B  but  with  glass  fiber  padding  substituted  for 
the  polyurethane  foam  of  the  seats.  Figure  3-41  shows  that  this 
substitution  resulted  in  significantly  less  heat  generation. 

The  agreement  between  the  computed  and  observed  temperatures 
is  good,  especially  near  the  end  of  the  test.  The  jagged 
appearance  of  the  computed  temperature  curve  reflects  the 
successive  ignition  and  burn  out  of  areas  of  material,  particularly 
the  seat  upholstery  fabric.  This  effect  is  noticeable  here 
because  of  the  absence  of  the  constantly  growing  seat  foam  fire 
that  was  present  in  the  earlier  cases.  Oxygen  consumption  was 
accurately  predicted  as  shown  in  Figure  3-42.  Agreement  in 
the  prediction  was  good  for  C02  but  the  concentration  of  the 
other  gases  was  overestimated  as  shown  in  Figures  3-43  through 
3-47.  Smoke  development  and  the  growth  of  the  upper  layer, 
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Figure  3-36. 
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Figure  3-37.  Hydrogen  Fluoride  Concentration  - Case  15B 
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TABLE  3-4 


COMPARISON  OF 
OF 

SIMULATED  AND 
DAMAGE  - CASE 

OBSERVED 

15B 

AREAS 

OBSERVED 

AREA111 

COMPUTED  AREA 

STRUCTURE 

(ft2) 

% TOTAL 
EXPOSED 

(ft2) 

% TOTAL 
EXPOSED 

Sidewall 

25.0 

29 

45.0 

100 

Ha track 

3.5 

5 

34.0 

100 

Ceiling 

59.0 

100 

86.0 

100 [2] 

Seats 

15.0 

33(3] 

60.5 

37[4] 

^ 1 ^ Damaged 
(center  7. 

area  within  test  section  considered 

5 ft  or  15  ft  length) 

by  model 

^ Difference  in  percent  figures  due  to  difference  in 
modeled  and  actual  surface  areas 

^ Damage  on  one  seat  row 

r 4 1 

1 'Total  damage  on  three  seat  rows 
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Figure  3-43.  Carbon  Dioxide  Concentration  - Case  15C 
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Figure  3-44.  Carbon  Monoxide  Concentration  - Case  15C 
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Figure  3-45.  Hydrogen  Chloride  Concentration  - Case  15C 
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Figure  3-46.  Hydrogen  Cyanide  Concentration  - Case  15C 


Figures  3-48  and  3-49,  resembled  the  earlier  cases.  Table 
3-5  and  Figure  3-50  show  the  predicted  and  observed  fire  spread. 

In  Figure  3-50  the  ceiling  and  hatrack  has  been  shaded  to 
show  fire  involvement,  although  in  the  test  only  heavy  sooting 
and  some  melting  or  burning  of  the  surface  coating  in  small, 
scattered  areas  was  the  only  evidence  of  fire  damage  on  these 
surfaces . 

3.3.6  Case  26P  - Present  In-Service  (1963)  Materials 

In  the  26-foot  fuselage  fire  tests,  much  more  tempera- 
ture and  gas  concentration  data  was  available  for  comparison 
to  the  model's  predictions.  Of  the  fifty  thermocouple  channels, 
four  were  selected  for  comparison  to  DACFIR's  computation  of 
the  upper  zone  gas  temperature  and  are  shown  in  Figure  3-51. 

The  broken  curves  represent  readings  at  12  inches  above  the 
floor  on  the  aft  cabin  tree,  48  inches  above  the  floor  on  the 
forward  cabin  tree,  48  inches  above  the  floor  on  the  center 
tree,  and  82  inches  above  the  floor  (8  inches  below  the  ceiling) 
on  the  center  cabin  tree.  All  readings  are  on  the  fuselage 
centerline.  The  solid  line  is  the  computed  temperature  which 
shows  a fair  agreement  with  the  48  inch  forward  thermocouple 
except  for  the  period  before  30  seconds.  The  analysis  of  this 
test  given  in  [6]  indicates  that  a flash  fire  occurred  around 
90  seconds  as  indicated  by  the  rapid  rise  of  the  ceiling  thermo- 
couple reading  at  that  time.  In  contrast  to  Case  15P,  however, 
the  temperature  data  shows  that  at  least  some  vigorous 
burning  continued  after  90  seconds.  Therefore  the  flash  fire 
was  probably  restricted  to  the  area  near  the  ignition  fire  and 
did  not  permeate  the  entire  cabin. 

Figures  3-52  through  3-58  give  the  cabin  gas 
composition.  Two  sets  of  experimental  data  are  included. 

Readings  at  sample  port  A were  taken  48  inches  above  the  floor 
and  directly  opposite  the  fire.  Readings  from  sample  port  B 
were  taken  at  the  same  vertical  position  and  side  of  the 
fuselage  as  A but  20  inches  from  the  forward  bulkhead.  Equip- 
ment failure  caused  much  of  the  gas  data  to  be  lost  after 
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Figure  3-49-  Smoke  Accumulation  - Case  15C 


TABLE  3-5 

COMPARISON  OF  SIMULATED  AND  OBSERVED  AREAS 
OF  DAMAGE  - CASE  15C 


OBSERVED 

AREA C 11 

COMPUTED 

AREA 

STRUCTURE 

(ft2) 

% TOTAL 
EXPOSED 

(ft2) 

% TOTAL 
EXPOSED 

Sidewall 

20.0 

23 

45.0 

100 

Hatrack 

0.0 

0 

33.8 

100 

Ceiling 

5.0 

4 

55.5 

64 

Seats 

19.0 

43<21 

124.2 

93[3 

Deunaged  area  within  test  section  considered  by  model 
(center  7.5  ft  or  15  ft  length) 

(21 

L 1 Daunage  on  one  seat  row 
^ Total  deunage  on  three  seat  rows 
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Figure  3-50.  Computed  and  Observed  Areas  of  Damage  - Case  15C 
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Figure  3-57.  Hydrogen  Fluoride  Concentration  - Case  26P 
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Figure  3-58.  Nitrogen  Dioxide  Concentration  - Case  26P 
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120  seconds.  For  0 2,  C02 , CO,  and  HF  the  agreement  of  the  model 
and  results  are  fair  to  good.  HCN  and  HC1  are  underestimated 
by  an  order  of  magnitude  at  some  points.  N02  levels  were  very 
low  in  the  test,  always  less  than  10  ppm,  and  the  computed 
value  of  N02  does  not  exceed  about  100  ppm. 

Figure  3-59  shows  that  the  computed  upper  zone  depth 
filled  the  2000  cubic  foot  volume  of  the  cabin  at  140  seconds. 
Figure  3-60  shows  that  the  smoke  reduced  visibility  below  1% 
at  115  seconds.  Table  3-6  and  Figure  3-61  show  the  comparison 
of  fire  involvement. 

3.3.7  Case  26N  - Improved  Materials 

Figures  3-62  through  3-71  and  Table  3-7  give  the 
comparison  between  the  model  and  test  results  for  Case  26N. 

In  this  test  improved  materials  were  used.  The  materials  were 
essentially  the  same  as  those  of  15B  but  with  the  addition  of 
carpet,  window  reveals  and  dust  panes,  and  simulated  PSU 
surfaces.  The  data  presented  is  that  from  the  second  attempt 
at  conducting  the  full  scale  test;  in  the  first  attempt  flames 
from  the  fuel  pan  would  not  enter  the  cabin  due,  apparently, 
to  wind  conditions  |6].  In  this  second  try,  the  same  problem 
with  ignition  occurred  until  the  fuel  pan  wind  screen  was 
removed  at  which  point  almost  all  the  flames  suddenly  entered 
the  rupture  driven  by  a six  knot  wind. 

The  larger  size  of  the  fuel  fire  and,  perhaps,  its 
intensification  by  the  wind  is  apparently  the  reason  that  the 
temperatures,  gas  concentrations,  and  areas  of  damage  in  this 
test  are  roughly  the  same  as  Case  26P  even  through  improved 
materials  were  used.  The  figures  do  not  show  the  improvement 
that  was  observed  between  Cases  15P  and  15B.  The  DACFIR 
model,  however,  demonstrates  the  improved  flammability  properties 
of  the  furnishings  as  is  shown  by  the  lower  computed  temperature, 
Figure  3-62  versus  Figure  3-51,  and  smaller  computed  flame 
spread.  Figure  3-71  versus  Figure  3-61.  Since  the  model  does 
not  contain  a mechanism  for  simulating  the  wind  effect  on  the 
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TABLE  3-6 


COMPARISON  OF  SIMULATED  AND  OBSERVED  AREAS 
OF  DAMAGE  - CASE  26P 


OBSERVED  AREA  ^ 1 ^ COMPUTED  AREA 


0 

% TOTAL 

% TOTAL 

STRUCTURE 

(ft:  ) 

EXPOSED 

(ft^) 

EXPOSED 

Carpet 

2.5 

3 

14.3 

17 

Sidewall 

25  - 34 

30  - 40[21 

43.3 

48 

PSU 

32.5 

100[2] 

14.3 

48 

Ha track 

10.7 

95 

2.8 

37 

Bullnose 

Ha track 

32.5 

100 

6.3 

5 

Ceiling 

75.0 

100  C3] 

13.8 

12 

Seats  . 

1 

0 

0 

0 

0 

2 

0 

0 

28.0 

51 

3 

0 

0 

0 

0 

4 

18.5 

42 

41.8 

76 

5 

0 

0 

0 

0 

6 

0 

0 

17.5 

32 

^ Damaged  area  within 

test  section 

considered  by 

model 

[2] 

Damage  described  as 

"melting  and 

shrinking"  of 

the  vinyl 

^ 3 ^ Much  delamination  occurred,  the  amount  of  actual  burning  is 
not  knovm 
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Figure  3-61.  Computed  and  Observed  Areas  of  Damage  - Case  26P 
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Figure  3-62.  Upper  Zone  Temperature  - Case  26N 
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Figure  3-63.  Oxygen  Concentration  - Case  26N 
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Figure  3-66.  Hydrogen  Chloride  Concentration  - Case  26N 
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Figure  3-67.  Hydrogen  Cyanide  Concentration  - Case  26N 
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Figure  3-68.  Hydrogen  Fluoride  Concentration  - Case  26N 
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Figure  3-69.  Lower  Zone  Depth  - Case  26N 
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TABLE  3-7 


COMPARISON  OF  SIMULATED  AND  OBSERVED  AREAS 
OF  DAMAGE  - CASE  26N 


OBSERVED  AREA^11  COMPUTED  AREA 


STRUCTURE 

(ft2) 

% TOTAL 
EXPOSED 

(ft2) 

% TOTAL 
EXPOSED 

Carpet 

2.0 

2 

Material 

Not  Available 

Sidewall 

25  - 34 

30  - 40 

13.5 

15 

PSU-Hatrack 

Bottom 

32.5 

100 

7.8 

52 

Hatrack 

Edge  & Top 

32.5 

100 

1.8 

2 

Ceiling 

75.0 [2] 

100 

5.8 

7 

Seats  SifHE 

0 

0 

0 

0 

2 

18.0 

40 

5.5 

10 

3 

0 

0 

0 

0 

4 

40.0 

90 

29.3 

53 

5 

0 

0 

0 

0 

6 

27.0 

60 

2.5 

5 

Damaged  area  within  test  section  considered  by  model 
(center  7.5  ft  or  15  ft  length) 


f 2 1 

Delamination  occurred,  the  amount  of  actual  burning  in 
not  known 
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fire,  the  comparison  between  the  computed  and  observed 
results  is  not  strictly  valid.  It  is  interesting  to  note,  however, 
that  the  agreement  for  gas  concentration  is  still  fair  during 
the  period  of  somewhat  steady  burning,  100  to  160  seconds  as 
indicated  by  the  temperature  readings.  For  this  larger 
cabin,  the  computed  upper  zone  did  not  reach  the  floor  by  the 
end  of  the  run.  Figure  3-69.  Smoke  generation  in  the  model 
was  quite  severe  as  shown  in  Figure  3-70. 

3.3.8  Summary  of  the  Model  and  Test  Comparisons 

The  agreement  between  the  results  of  the  seven  AIA- 
CDP  cabin  mock-up  fire  tests  and  the  predictions  of  the  DACFIR 
model  are  summarized  in  Table  3-8.  The  table  shows  that  the 
program  was  most  successful  at  predicting  the  order  of  magnitude 
of  the  accumulation  of  toxic  gases.  Temperature  prediction  was 
generally  good  for  later  times  in  the  tests  but  never  very 
good  for  the  early  times.  Prediction  of  the  rate  of  fire  spread 
and  material  involvement  was  usually  fair  or  good  for  the  overhead 
surfaces,  the  ceiling,  hatrack,  and  upper  sidewall,  but  poor  for 
the  lower  sidewall,  floor,  and  seats.  Overestimation  of  the 
spread  of  flames  over  these  surfaces  could  be  due  to  the  flame 
radiation  model  and/or  problems  with  the  flame  spread  rate 
input  data.  Since  the  production  of  heat,  smoke,  and  toxic  gas 
is  a function  of  the  area  of  each  material  burning,  improvement 


in  the  spread  prediction  should  result  in  improvement  in  the 
temperature  and  gas  concentrations  computed. 


77 


SUMMARY  OF  MODEL  AND  TEST  COMPARISONS 


■ 
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Third  100  seconds  of  test 


SECTION  4 
CONCLUSIONS 

After  consideration  of  the  comparisons  of  the  DACFIR 
model  to  the  AIA-CDP  cabin  fire  tests,  the  following  conclusions 
have  been  reached  concerning  the  performance  of  the  model  and 
the  nature  of  the  input  data. 

4.1  SIMULATION  OF  SPECIFIC  TESTS 

The  DACFIR  model's  ability  to  reproduce  the  results  of  the 
AIA  tests  has  been  summarized  in  Table  3-8.  From  examination 
of  the  table,  it  can  be  seen  that  the  model  has  been  moderately 
successful  in  predicting  the  average  temperature  in  the  upper 
portion  of  the  cabin.  Better  results  occur  for  the  smaller 
cabin  size,  as  might  be  expected  by  the  assumption  of  a 
uniform  hot  gas  zone.  The  model  usually  does  best  at  predicting 
toxic  gas  accumulation,  in  the  sense  of  an  order  of  magnitude  at 
least.  The  model  does  not  seem  to  do  well  at  predicting  damage 
(areas  of  fire  spread  or  other  thermal  damage)  except  for  the 
overhead  surfaces.  Overall,  the  program  can  distinguish  the 
effects  of  relatively  more  flammable  or  less  flammable  furnishings. 

Specific  findings  and  conclusions  about  the  assumptions 
and  structure  of  DACFIR  and  about  the  laboratory  data  used  as 
its  input  are  listed  below. 

1.  The  assumption  of  a uniform  temperature  zone 
of  combustion  products  (the  upper  zone)  is  probably  an  over- 
simplification. Significant  horizontal  and  vertical  temperature 
gradients  can  occur  in  the  cabin  as  is  apparent  in  Cases  26P  and 
26N.  The  model  should  be  refined  to  account  for  these  variations 
to  better  model  heat  transfer  to  objects  remote  from  the  fire  as 
well  as  for  the  effects  on  passenger  evacuation.  It  should  be 
realized  that  some  thermocouple  readings  from  the  burn  tests  in- 
dicate that  the  thermocouple  was  within  or  very  close  to  the 
flames  of  the  fuel  and  seat  fires  and  should  not  be  compared  to  the 


average  upper  zone  temperature  computed  by  DACFIR.  Neverthe- 
less, a better  gas  flow  and  thermodynamics  model  should  be 
sought.  A promising  improved  model  is  presented  in  Appendix  A. 


2.  The  burning  rate  of  the  seat  material  as 
observed  in  the  full  scale  tests  does  not  appear  to  be  constant 
over  the  life  of  a burning  element  as  is  assumed  by  the  program. 
The  burning  rate  seems  to  accelerate  significantly  leading 

to  more  rapid  release  of  unwanted  quantities  and  depletion  of 
oxygen.  Modification  of  the  DACFIR  program  to  remove  the 
assumption  of  constant  release  rates  of  smoke,  heat,  and 
gases  (at  constnat  applied  flux  levels)  is  possible.  The 
modifications  necessary  are  straightforward  but  would  increase 
significantly  the  volume  of  input  data  required.  Fortunately, 
data  in  the  required  form  is  already  available  from  past 
laboratory  test  programs. 

3.  With  regard  to  the  laboratory  test  data,  it 
was  determined  that  there  appears  to  be  a significant  effect  in 
the  heat  release  rate  measurements  of  the  thermal  inertial  of 
the  OSU  analyzer.  This  was  indicated  by  the  disagreement  of  the 
time  to  ignite  and  time  to  burnout  measured  visually  and  from 
the  heat  release  curve.  If  the  thermal  inertia  effects  can 

be  removed  from  the  heat  release  data,  the  result  could  be  more 
vigorous  heat  release  rates  which  would  make  the  model  show 
more  rapid  burning  development,  as  seems  to  be  the  case  with  the 
seat  materials  in  the  tests. 

4.  The  DACFIR  model  consistently  over  predicts 
the  rate  and  extent  of  flame  spread,  particularly  on  vertical 
surfaces.  The  problem  may  be  due  to  one  or  both  of  the  follow- 
ing factors.  The  method  of  measurement  of  the  flame  spread 
velocity  as  a function  of  applied  heat  flux  in  the  OSU  analyzer 
may  not  be  reliable,  or  the  radiation  feedback  intensity 
calculation  used  in  the  model  to  select  the  flame  spread 
velocity  may  be  over  estimating  the  radiant  level  at  the  edge 
of  the  flame.  In  addition,  the  representation  in  the  model  of 
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[seat  rows  as  a continuous  bench-like  arrangement  allows  fires 
to  spread  over  an  entire  row  without  hindrance,  whereas  in 
the  tests  the  breaks  between  individual  seats  seem  to  stop  or 

I slow  the  fire  spread. 

5.  The  behavior  of  the  ignition  source  requires 
further  study.  The  importance  of  the  ignition  source  is 
demonstrated  by  the  outcome  of  Case  26N  where  the  ignition 
source  size  and  the  effect  of  the  wind  overcame  whatever 
improvement  was  to  be  obtained  by  the  new  materials.  Proper 
understanding  of  the  furnishing  materials'  effect  on  the  cabin 
environment  during  a fire  requires  a better  understanding  of 
the  character  of  the  ignition  source. 

6.  For  the  more  vigorous  fires,  oxygen  depletion 
affects  fire  behavior  in  the  later  stages  of  the  test.  While 


the  DACFIR  model  computes  oxygen  depletion,  the  mechanism 
coupling  the  lower  level  of  oxygen  to  fire  behavior  is  not 
included  due  to  lack  of  an  available  quantitative  relationship. 

4.2  UNDERSTANDING  OF  THE  TESTING/MODELING  PROCESS 

The  development,  comparison,  and  refirment  of  the  DACFIR 
model  has  proven  to  be  a great  aid  in  organizing  and  evaluating 
the  results  from  full-scale  tests.  By  attempting  to  model  a 
full-scale  test  mathematically,  the  fire  and  its  development 
must  be  viewed  as  a system  of  coupled,  interdependent  processes. 
Each  process  - the  flame,  the  materials'  thermal  degradation, 
the  gas  dynamics,  etc.  - must  be  understood  to  the  precision 
necessary  for  the  whole  mathematical  system  to  reproduce  the 
full-scale  results.  In  building  and  testing  the  complete 
model,  the  relative  importance  of  each  part  can  be  found.  This 
knowledge  can  then  be  applied  to  the  design  of  new  tests  to 
determine  which  situations  (scenarios)  are  important,  what 
processes  must  be  carefully  controlled,  and  what  measurements 
and  accuracies  are  required. 


APPENDIX  A 

AN  IMPROVED  CABIN  GAS  DYNAMICS  MODEL 

INTRODUCTION 

The  comparison  of  the  DACFIR  model  to  the  AIA-CDP  fire  tests 
indicates  that  the  simple,  two  zone  treatment  of  interior  gas 
dynamics  may  not  predict  well  the  observed  temperature  field  in 
cabin  fires  for  certain  situations.  It  is  apparent  that  significant 
vertical  and  horizontal  temperature  gradients  can  exist  within 
the  region  of  smoke  and  gas  accumulation  if  the  cabin  is  much 
longer  than  its  width.  Since  actual  cabins  have  a greater 
length  to  width  ratio  than  the  mock-up  cabins  used  in  the  com- 
parisons, a method  is  needed  to  improve  the  model's  gas  dynamics. 

Few  researchers  have  dealt  directly  with  the  flow  of  hot  gas 

under  a horizontal  ceiling.  Among  those  that  have  are  Alpert  ^ 

[21 

and  Hwang,  et  al.  Alpert  analyzed  the  ceiling  jet  resulting 

from  the  weakly  buoyant  fire  plume  striking  a flat  ceiling. 

His  model  assumes  an  infinite  ceiling  extent  and  uses  an 
axisymmetric  geometry.  Hwang  and  his  co-workers  studied  the 
phenomenon  of  reverse  flow,  i.e.,  jet  flow  against  a 


ventilating  current,  in  mine  shaft  fires.  They  formulate  a 
quasi-steady,  two-dimensional  model  applicable  to  fire  gas 
flows  in  tunnel-like  geometries.  Naturally,  this  is  suggestive 
of  the  aircraft  problem  and  has  motivated  the  work  detailed  below. 

The  feasibility  of  applying  a horizontal  jet  flow  model  to 
aircraft  fires  was  investigated  in  the  simplified  case  of  a 
single  plume  generating  a jet  unconfined  by  end  walls. 

Theoretical  predictions  have  been  compared  to  results  from  a 
series  of  cabin  fire  tests  conducted  in  a B737  fuselage  at 
Johnson  Space  Center  of  the  National  Aeronautics  and  Space 


1 JAlpert,  R.L.,  "Turbulent  Ceiling-Jet  Induced  by  Large-Scale  Fires", 
Combustion  Science  and  Technology,  1975,  Vol.  11,  pp.  197-213. 

1 1 Hwang,  C.C.  et  al.,  "Reverse  Stratified  Flow  in  Duct  Fires:  A 
Two-Dimensional  Approach",  Sixteenth  Symposium  (International)  on 
Combustion  (1976),  p.  1385. 
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Administration  ( NASA- LB J) . A description  of  the  model  and  the 
results  of  this  first  test  are  given  below. 

THE  CEILING  JET 

Figure  A.l  presents  the  features  of  the  jet  model.  The 
vertically  rising  plume  from  a fire  near  the  cabin  floor  strikes 
the  flat  and  level  ceiling.  Plumes  in  the  DACFIR  model  are 
axisymmetric . The  problems  of  determining  the  flow  patterns  in 
the  turning  region  - where  the  plume  gas  turns  and  eventually 
assumes  a more  or  less  parallel  flow  away  from  the  fire  - are 
not  thought  worth  the  considerable  effort  for  solution  at  present. 
Therefore,  we  assume  that  the  plume  flow  into  the  turning 
region  is  immediately  redirected  and  redistributed  into  a starting 
flow  for  the  ceiling  jet  uniform  across  the  cabin  width.  In 
this  model  the  cabin  is  assumed  to  be  rectangular  in  cross 
section.  The  figure  shows  the  assumed  velocity  profile 
in  the  ceiling  jet.  The  maximum  gas  velocity  and  minimum  density 
(and  thus  the  maximum  temperature)  occur  at  the  ceiling.  The 
decrease  in  velocity  and  increase  in  density  through  the  jet  is 
taken  to  be  semi-gaussian  with  the  lower  edge  of  the  jet 
selected  as  the  point  where  the  velocity  has  fallen  to  five 
percent  of  its  maximum  value.  The  jet  exchanges  momentum  and 
heat  with  the  ceiling  by  shear  and  convective  heat  transfer. 

At  the  lower  edge,  turbulent  entrainment  of  the  cooler  air  on 
which  the  jet  floats  carries  in  mass,  momentum,  and  energy. 

The  cool  lower  air  may  be  flowing  toward  the  fire  driven  by 
the  entrainment  of  the  plume  and  also,  if  applicable,  by  a 
ventilation  system. 

Following  Hwang,  et  al.,  the  integral  formulation  for 

<» 

the  equations  of  motion  of  the  ceiling  jet  are: 

-A-  ovdy  = aDa,(V  + V„> 

’I  <lvCPWy  * cpT.».“<'?  + vJ-%, 


(1)  Continuity 

(2)  Energy 


A-2 


(3)  Momentum  (in  the 

horizontal  direction) 


d 

dx 


/oh  pv2dy 


-op  V (V  + VJ  - 


d 

dT 


yg(pa>-p)dy-Tw-Te./oh(|Z)a)dy 


The  coordinate  system  of  this  model  has  its  origin  at  the 
ceiling  impact  point  with  X the  horizontal  coordinate  (positive 
X directed  against  the  ambient  flow)  and  Y the  vertical  coordinate 
(positive  Y directed  downward) . 

The  variables  are  density  (p)  , velocity  (v) , and  thickness 
(h)  of  the  ceiling  jet  layer.  Other  quantities  are: 

a entrainment  parameter,  a constant  for  this  model; 

C_  specific  heat  of  lower  air  (ventilation  flow)  and  of 
the  combustion  products  (assumed  constant) ; 

(-2— ) pressure  drop  in  ventilation  flow; 

□ X 00 

gravity; 

convection  heat  loss  to  the  ceiling; 
shear  stress  at  layer  interface; 
wall  shear  stress  at  the  ceiling; 
temperature;  and 

characteristic  velocity  of  the  jet. 

The  subscript  (°°)  denotes  ambient  conditions. 

Once  profiles  for  density  and  velocity  are  specified,  the 
integrals  in  equations  (1)  - (3)  may  be  evaluated  and  the  resulting 
set  of  ordinary  differential  equations  solved. 

The  semi-guassian  profiles  assumed  are: 

v = vmexp(-y2/i2)  p * pw-Ap  = P^-Ap^xp (-y2/?.2) 

where  l is  the  profile  half  width.  To  relate  l to  the  jet 
thickness  h,  the  thickness  is  defined  to  be  at  the  point  on  the 
velocity  profile  where  the  velocity  has  fallen  to  five  percent 
of  its  maximum  value.  This  gives 

l = h//3 
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With  these  substitutions,  equations  (1)  - (3)  may  be 
restructured  into  the  following  set: 
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where  T is  the  wall  temperature,  h^  is  the  convection  coefficient 
w t 

(assumed  constant)  , T is  the  jet  temperature  (=P00Too/p)  , and  Cf  is 
a friction  factor,  taken  as  = 0.008. 

Thus  the  set  is  linear  in  the  derivatives  and  may  be 
integrated  easily. 

The  solution  produces  a longitudinal  density  gradient;  hence 
the  desired  temperature  gradient.  Since  a density  profile 
is  assumed,  a vertical  gradient  is  also  calculable.  Consequently, 
it  is  possible  to  impose  a grid  on  the  cabin  interior  and  thereby 
obtain  a set  of  temperatures  from  which  the  temperature  at  a given 
point  may  be  interpolated. 

In  addition,  application  of  conservation  of  species  allows 
the  calculation  of  smoke  and  toxic  gas  concentration  gradients  : 


(4)  _d_  -h  pvY.dy  = aP<x)  (V  + VJ  Y^. 

dx  o 

where  Y.  is  the  mass  fraction  of  ith  gas  or  smoke  in  jet  layer, 

i 4*  j>\ 

and  Y . is  the  mass  fraction  of  i gas  or  smoke  in  ambient  layer. 
00  1 

If  uniform  mixing  within  the  hot  layer  is  assumed,  along  with 
the  previous  density  and  velocity  profiles,  equation  (4)  becomes 
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INITIALIZATION  BY  RISING  FIRE  PLUMES 

To  initialize  the  ceiling  jet,  the  Steward-Fang  plume  is  used 
to  calculate  upward  mass  and  heat  flows  from  a fire  in  the  lower 
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cabin.  If  in  the  turning  region  (where  upward  flow  impinges  the 
ceiling  and  is  turned  horizontal)  heat  and  momentum  loss  are 
assumed  negligible,  then  the  conservation  laws  imply  the  following 
simple  relations  between  plume  and  initial  jet  values. 


(5) 

Mass 

2 h ' 

Trr  p v = 2W  / p .v,  dy 

P P P o j 

(6) 

Momentum 

Trr  p v = 2W  /n  p .v;  dy 

P P P o j 1 1 

(7) 

Energy 

irr  2C  T p v = 2W  Z*1  C T p v.dy 

P PP'PP  O P oo^oo  j ■* 

(8) 

z - h ' = y 
•'o 

where 

z is  the  cabin 

height,  h'  is  the  initial  jet  thickness, 

is  the  distance  to  the  jet  layer  from  the  floor,  and  W is  the 
width  of  the  turning  region.  Subscripts  p and  j indicate  plume 
and  jet  values  respectively.  The  factor  of  two  occurs  in  the 
right-hand  side  of  equations  (5) , (6) , and  (7)  since  it  is  assumed 

that  the  upward  flow  divides  equally  between  the  positive  and  nega- 
tive X directions. 

With  the  assumed  jet  profiles  of  density  and  velocity,  equa- 
tions (5)  - (7)  become: 


(5a)  ir  2p  v = W /tT(o  -io  ) erf  (/3)//3  h'v 

p p p °°  m m 


With  the  model  in  this  form,  comparisons  were  made  with  one  of 
the  NASA-LBJ  "Design  Fire"  tests  in  the  737  fuselage^.  Starting 
with  an  assumed  mass  loss  and  heat  release  rate  for  the  24  X 24 
inch  square  pan  of  Jet-A  used  in  the  test,  computations  were  made 
for  the  hot  gas  flow  from  the  ceiling  impact  point  to  a position 
30  feet  down  the  fuselage.  A ventilation  flow  velocity  of 
8 ft/min  was  used  which  corresponds  to  the  475  cfm  ventilation 
rate  used  in  the  burn  test. 

Initial  results  of  the  comparison  indicated  that  the  entrain- 
ment constant  for  the  ceiling  jet  was  too  small  to  account  for 
the  temperature  decrease  with  distance  observed  at  NASA-LBJ. 

Some  adjustment  of  the  constants  in  the  Steward-Fang  plume  model 
were  also  necessary  to  bring  the  temperature  at  the  plume  impact 
point  up  to  the  observed  value.  After  these  adjustments,  the 
performance  of  the  model  is  shown  in  Figure  A. 2. 

In  the  figure,  the  computed  gas  temperatures  at  the  ceiling  and 
four  inches  below  the  ceiling  are  shown  by  the  solid  and  broken 
curves  respectively.  Two  curves  of  each  type  are  shown,  one  for 
the  flow  against  the  ventilation  current  ("upwind")  and  one  for 
flow  in  the  same  direction  as  the  ventilation  current  ("downwind"). 
The  square  and  circular  symbols  are  thermocouple  readings  from 
NASA-LBJ  test  number  17.  The  square  symbols  are  values  downwind 
from  the  flow  and  the  circles  are  upwind  readings.  All  thermo- 
couples were  four  inches  below  the  ceiling. 

The  figure  shows  that  the  computed  temperature  at  the 
ceiling  agrees  very  well  with  the  measured  temperature  four  inches 
below  the  ceiling.  Except  for  one  point,  the  computed 
temperature  four  inches  below  the  ceiling  is  consistently  lower 
than  the  measurements.  Despite  this  discrepancy  in  the  vertical 
direction,  the  measured  and  computed  horizontal  temperature  changes 
are  in  excellent  agreement.  Although  only  one  case  was  analyzed 
and  some  room  for  improvement  has  been  noted,  the  results  shown 
here  are  encouraging. 


(3]  Tustin,  E.  A.,  et  al,  "Development  of  Fire  Test  Methods 
for  Airplane  Interior  Materials,"  Boeing  Commerical  Airplane 
Co.,  NASA-CR,  in  preparation. 
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Distance  from  Fire  (ft) 

Figure  A-2.  Comparison  of  Horizontal  Gas  Flow  Model  to  NASA 
737  Fire  Test  17  (2-1-77) 


APPENDIX  B 

LABORATORY  DATA  COLLECTION 

This  appendix  consists  of  the  final  technical  report  on 
the  laboratory  test  data  collection  program  conducted  in  support 
of  the  DACFIR  validation  effort. 

A very  large  amount  of  data  was  collected  consisting  of 
flame  spread  rates;  time  constants;  smoke,  heat,  and  gas  release 
rates;  and  weight  loss  measurements.  Since  many  of  the  cabin 
materials  tested  are  either  no  longer  in  widespread  use  or  were 
experimental  at  the  time  of  the  full  scale  tests,  it  was  not 
deemed  justified  to  present  the  data  in  this  report.  The  data 
can  be  made  available  to  interested  parties  by  contacting  the 
authors  or  the  FAA  technical  monitor. 
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Testing  by  The  Boeing  Company  in  response  to  the  University  of  Dayton 
P.0.  #RI-77086  dated  7-9-76  for  "Verification  of  the  DACFIR  "Model ing 
Concepts  Using  Data  from  Past  Burn  Tests"  has  been  completed.  Incremental 
shipments  of  the  total  information,  photographs,  16-mm  color  films,  and 
data  have  been  made  from  time  to  time  in  care  of  Mr.  Charles  MacArthur, 
University  of  Dayton. 


DACFIR  LABORATORY  TEST  PROGRAM 


SUMMARY 

The  testing  reported  herein  is  in  response  to  the  work  requested  in 
the  University  of  Dayton  Research  Institute  P.0.  #RI-77086  dated  7-9-76. 

Location  of  some  of  the  materials  used  in  the  1968  fire  testing  pro- 
gram was  responsible  for  a large  part  of  the  delay  in  getting  the  labora- 
tory work  completed.  Many  of  these  materials  are  no  longer  used,  some  were 
never  used,  and  others  have  become  obsolete  with  material  suppliers, 
i Therefore,  while  some  materials  were  not  available  for  testing,  a best 
effort  possible  on  the  part  of  Boeing  has  been  followed  in  obtaining  data 
to  satisfy  the  intent  of  the  program. 

DESCRIPTION  OF  MATERIALS 

All  of  the  31  materials  requested  in  the  University  of  Dayton  Research 
Institute's  letter  dated  10  August  1976  were  accumulated  and  tested  in  the 
National  Bureau  of  Standards  (NBS)  and  Ohio  State  University  (OSU)  test 
chambers,  except  the  following: 

Material  No.  3 --  Enameled  aluminum  with  Tedlar  surface  film 

Material  No.'  24  --  Muslin-covered  polyurethane,  upholstered  in  cotton 

fabric 

Material  No.  28  --  Olive  drab  Nomex 
These  materials  could  not  be  located  or  reasonably  simulated. 
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However,  it  should  be  noted  that  duplicate  materials  were  not  tested. 
Therefore,  there  are  no  data  transmitted  for  the  following: 

Material  No.  5 and  Material  No.  12  (same  as  No.  4) 

Material  No.  10  (same  as  No.  8) 

The  materials  tested,  description,  material  use  (tests),  and  specimen 
construction  were  as  indicated  in  table  I. 

OSU  COMBUSTION  ANALYZER  TESTING 


All  materials  were  tested  at  four  heat  flux  levels  (5.0,  3.5,  2.5,  and 
1.5  watts/sq  cm),  with  orientation(s)  as  indicated  in  the  University  of 
Dayton  letter.  Due  to  the  limitations  of  materials,  time,  or  anticipated 
results,  only  one  specimen  was  tested  at  the  lower  two  heat  flux  levels. 
Three  specimens  per  material  were  tested  at  the  two  higher  heat  flux  levels. 


The 

following  is  the  test  procedure 

used  for  the  OSU  testing: 

0 

Specimen  weighed 

Recorded 

0 

Specimen  inserted  into  chamber 

Time  noted 

0 

Start  of  smoke  (visual) 

Time  measured  (seconds) 

0 

Start  of  flame  (visual) 

Time  measured  (seconds) 

0 

Spread  of  flame  (visual) 

Time  measured  (seconds) 

0 

End  of  flame  (visual) 

Time  measured  (seconds) 

0 

End  of  smoke  (visual ) 

Time  measured  (seconds) 

0 

Specimen  removed  from  chamber 

Time  noted 

0 

Specimen  weighed 

Recorded 

0 

Strip  chart 

Recorded  heat  and  smoke 

curves 

The  smoke  and  heat  curves  were  digitized  with  a direct  link  to  a CDC  6600 
computer. 
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CALCULATIONS  AND  MEASUREMENTS 


The  following  calculations  and/or  measurements  were  obtained  from 
tests  made  in  the  OSU  combustion  analyzer: 

o Time  to  flame  (t^)  HRC  --  Time  between  specimen  insertion  and  the 
time  the  heat  release  rate  curve  begins  to  rise  above  the  zero 
1 evel . 

o Time  to  flame  (t^)  visual  — Time  between  specimen  insertion  and 
visual  appearance  of  flame. 

o Flame  spread  rate  — Time  required  for  flames  to  travel  from  the 
point  of  flame  application  to  a specific  edge  of  the  specimen. 

In  the  case  of  horizontal  test  position  for  flaming  condition,  a 
single  point  of  flame  was  always  applied  at  the  center  of  the 
specimen  and  the  rate  of  spread  was  taken  over  the  5-inch  dis- 
tance to  the  farthest  edge.  For  vertical  test,  a multiple  flame 
was  applied  at  edge  or  top  or  bottom  and  the  rate  of  spread  was 
taken  over  the  6-inch  distance  to  the  opposite  edge, 
o Time  to  char  (t^c)  visual  --  Time  between  first  appearance  of 
flame  and  flame  out  or  ceased  smoking, 
o Time  to  char  (tf(.)  HRC  — Unobtainable.  (In  the  interest  of 
reducing  overall  test  time,  the  chamber  was  not  permitted  to 
return  to  the  zero  level  before  removal  of  specimen  from  chamber.) 
o Total  heat  released  (ht)  (Btu)  — Computed  from  the  heat  release 
rate  curve  generated  by  the  OSU  combustion  analyzer. 


0 


Total  smoke  released  in  the  flaming  state  (DSf)  and  in  the  smol- 
dering state  (Ds$)  — Computed  from  the  smoke  release  curves 
generated  by  the  OSU  combustion  analyzer.  Also,  a determination 
was  made  of  the  maximum  value  of  the  specific  optical  density 
from  the  smoke  curve. 

Weiaht  loss  rate  = weight  before  - weight  after  . grams., 
weight  loss  rate  time  tQ  burn  minute 


SIZE  OF  SPECIMENS  FOR  OSU  TESTS 


o Horizontal  tests  --  4 inches  x 10  inches 

o Vertical  tests  — 6 inches  x 6 inches  (except  in  the  case  of 

Material  No.  14,  in  which  the  specimen  size  was  3 inches  x 3 
inches) . 

TOXIC  GAS  EVOLUTION  TESTING 


The  materials  listed  in  table  I were  burned  or  pyrolyzed  in  an  NBS 
Smoke  Chamber.  The  specimen  configuration  (nominal  3-inch  square),  com- 
bustion chamber  volume  (516  liters),  pilot  flame  size,  chamber  temperature, 
and  sample  exposure  conditions  were  all  as  described  for  the  standard  NBS 
Smoke  Density  test  (ref.  1:  NBS  Technical  Mote  708,  "Interlaboratory 
Evaluation  of  Smoke  Density  Chamber,"  T.  G.  Lee,  National  Bureau  of  Stan- 
dards, 1971).  Smoke  densities  were  measured  for  these  materials  but  are 
not  reported  under  this  contract. 
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The  evolution  of  certain  toxic  gases  (HC1,  HF,  HCN , riOx , SOg)  was 
measured  using  gas  detector  tubes  located  in  the  center  of  the  chamber 
above  the  specimen.  At  specific  times  during  the  test  (1,  2,  4,  and  10 
minutes),  aliquot  portions  of  the  chamber  atmosphere  were  drawn  through  the 
detector  tubes,  using  a semi-automatic  pumping  system  located  outside  of 
the  chamber.  Following  the  test,  the  tubes  were  removed,  and  the  concen- 
trations of  the  specific  toxicants  were  calculated  from  the  length  of  the 
observed  color  change  (stain)  produced. 

The  evolution  of  two  gases  was  monitored  by  pumping  chamber  gas 
through  two  specific  analytical  instruments,  a Beckman  Model  865  Nondisper- 
sive  Infrared  Analyzer  for  CO,  and  a Wilks  Instrument  Co.  Mi  ran  I Infrared 
Analyzer  for  CO2.  Concentrations  of  these  gases  at  1,  2,  4,  and  10  minutes 
were  read  from  the  continuous  recorder  traces. 

Specific  gas  detector  tubes  manufactured  for  Dragerwerks  A.G.  of 
Lubeck,  Germany  (Drager  tubes)  were  used  for  all  gases  except  HC1 , which 
was  determined  using  a detector  tube  manufactured  by  Mine  Safety  Appli- 
ances. Tests  were  conducted  only  for  toxicants  whose  formation  was  anti- 
cipated, based  upon  the  chemical  composition  of  the  test  material. 

Gas  detector  tubes  were  developed  for  determining  specific  toxicants 
in  clean  atmospheres,  i.e.,  as  in  industrial  hygiene  applications.  Accur- 
acies of  about  + 20  percent  are  typically  claimed  by  the  manufacturer. 
Although  they  have  been  widely  applied  for  analysis  of  combustion  gases  in 
the  NBS  Smoke  Chamber  and  elsewhere,  this  application  is  subject  to  the 
possibility  of  unknown  interferences.  In  our  opinion, vthe  results  obtained 
by  this  method  are  probably  accurate  within  + 50  percent. 
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1 


Gas  evolution  measurements  were  made  at  flux  levels  of  5.0,  3.5,  2.5, 
2 , 


and  1.5  watts/cm  in  both  the  flaming  mode  (sample  impinged  by  a pilot 
flame)  and  in  the  smoldering  mode  (sample  irradiated  without  an  ignition 
source  present).  Cost  and  time  limitations  prohibited  duplicate  determin- 


ations of  each  toxicant  at  1,  2,  4,  and  10  minutes  at  each  flux  level. 

.2 


Testing  was  begun  at  5.0  watts/cm  and  proceeded  to  lower  flux  levels.  As 
soon  as  the  amount  of  a specific  toxicant  dropped  to  a negligible  value, 
testing  for  that  toxicant  was  discontinued,  and  when  the  evolution  of  all 
toxicants  had  dropped  to  negligible  values,  testing  for  that  material  was 
discontinued.  For  this  reason,  very  few  data  were  collected  at  1.5  watts/. 
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APPENDIX  C 
DACFIR2  USER'S  GUIDE 


This  appendix  is  a guide  for  use  of  Version  2 of  the 
Dayton  Aircraft  Cabin  Fire  simulation  program  (DACFIR2). 

DACFIR2  differs  substantially  from  Version  1 of  the  DACFIR 
program,  a user's  guide  for  which  was  given  in  Volume  III 
of  Reference  [1] . While  this  appendix  draws  from  material 
presented  in  the  earlier  volume,  familiarity  with  Version  1 
of  DACFIR  is  not  assumbed  for  the  use  of  the  present  material. 

The  intent  of  this  user's  guide  is  to  provide  instruction  for 
the  efficient  use  of  DACFIR2,  but  not  to  present  the  construction 
of  the  computer  code  in  detail.  The  guide  provides  an  annotated 
flow  chart  of  the  main  controlling  program,  instructions  for 
preparing  the  input  data,  sample  input  and  output,  program 
statistical  data,  and  information  for  obtaining  copies  of  the 
program.  Due  to  the  large  size  of  DACFIR,  over  three  thousand 
source  statements,  a listing  of  the  code  is  not  included. 

C.l  BASIC  DEFINITIONS  AND  CONVENTIONS 

In  the  DACFIR  program,  cabin  interior  geometry  is  represented 
by  a connected  group  of  horizontal  and  vertical  surfaces.  The 
program  recognizes  a maximum  of  twenty  (20)  cabin  lining 
surfaces  (floor,  sidewalls,  ceiling,  etc.)  and  nine  (9)  seat 
groups.  The  cabin  lining  surfaces  and  seat  surfaces  are 
divided  into  square  elements  of  fixed  dimension,  0.5  by  0.5  ft, 
for  purposes  of  tracking  the  fire  growth.  Each  surface  and 
seat  group  is  identified  by  a single  number  and  each  element  on 
a surface  by  a pair  of  numbers,  the  element  i,j  indices. 

Assignment  of  the  surface  and  seat  numbers  and  element  indices 
is  made  by  the  program  based  upon  the  values  of  the  input  data 


^ Reeves,  J.B.,  and  C.D.  MacArthur,  "Dayton  Aircraft  Cabin  Fire 
Model,"  Volumes  I,  II,  and  III,  FAA-RD-76-120 , June  1976.  Volume 
III,  "Computer  Program  User's  Guide"  was  written  by  P.M.  Kahut. 
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describing  the  surfaces  and  seat  groups  and  upon  the  order  of 
the  input  data.  DACFIR2  is  designed  to  minimize  the  amount  of 
geometric  input  information  required.  To  do  so,  a number  of 
conventions  and  assumptions  are  adopted  of  which  the  user  must 
be  aware  in  preparing  the  geometric  input.  The  conventions  and 
assumptions  are  as  follows. 

(1)  All  cabin  lining  surfaces  are  assumed  to 

be  parallel  to  the  cabin  y-axis  (See  Figure 
C-l  and  Figure  2-1,  Section  2)  and  each 
surface  must  be  eitherparallel  to  perpen- 
dicular to  the  x-y  plane  (floor  plane) . 

The  seat  group  configuration  is  as 
described  in  Section  2 , all  dimensions 
fixed  except  width,  and  all  seats  must 
face  forward. 

(2)  All  cabin  lining  surfaces  are  assumed  to 
extend  unbroken  from  the  start  of  the 
"detailed"  section  length  (Figure  2—1 ) 
to  its  end;  this  distance  is  fixed  as 
7.5  ft,  that  is  15  element  lengths. 

(3)  Numbering  of  the  cabin  lining  surfaces 
always  starts  with  the  floor  as  surface 
number  one  and  proceeds  counterclockwise 
(facing  aft)  across  the  floor,  up  the  side- 
wall,  across  the  ceiling  and  down  the 
sidewall  to  return  to  the  floor.  This 
scheme  is  shown  in  Figure  C-l. 

(4)  Each  seat  group  is  constructed  of 
seven  surfaces  shown  in  Figure  C-2. 

Numbering  of  the  elements  of  a seat 
group  is  fixed  by  the  program  and  starts 
with  the  cushion  bottom,  proceeding  up 
the  backrest,  over  the  backrest  top,  down 
the  front  to  the  seat  cushion  top,  over 
the  cushion  top  and  returning  to  the  front 
edge  of  the  cushion  bottom. 
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i=23 


Figure  C-l.  Cabin  Lining  Surfaces.  Configuration  shown 
is  that  used  for  Cases  26P  and  26N. 
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Figure  C-2.  Seat  Surfaces 


C.2  PROGRAM  CONSTRUCTION 


The  following  pages  present  an  annotated  flow  chart  of 
the  main  controlling  program  of  DACFIR2 . The  flow  chart 
shows  the  order  of  call  of  each  subroutine  designed  for  a 
specific  task  or  process.  Comments  on  the  flow  chart  are 
identified  by  the  numerals  in  the  left  margin.  The  list  of 
comments  follows  the  flow  chart. 
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'MAIN  PROGRAM'  FLOW  CHART 


Comment  No. 


1. 


2. 

3. 

4. 

5. 

6. 

7. 

8. 

9.  0 


10. 


START 

1 

NPASS  = 0 
IPRl  - 0 
IPR2  = 2 

I 

CALL  INPUTG 
(IN) 


(IN  = 0)? — PRINT ; 

yes 

I 

PRINT:  'RUN  ID' 

I 

CALL  INPUTM 

CALL  IN IT 
CALL  INIT2 

I 

CALL  INPUTO 

1 

CALL  PRTINP 

I 

ITIM2  - ITIM2  + ITSPRD 
1 = 0 
K - 0 
KPR  = 0 
ISW  - 0 
ISAVE  = 0 

ITIME  - ITIME  - IDELT 

I 

TDQ  - 0. 

TRSF  « 0. 

TOTVIT  - 0. 

TRGF(J)  - 0;  J - 1,10 


'NORMAL 
COMPLETION ' 


STOP 


'MAIN  PROGRAM'  FLOW  CHART  (Continued) 


Comment  No. 


1 = 1 + 1 


NFE(I)  = 0? 


■CALL  FIRE (I) 

I 

K = KPR?  — 


CALL  RATES (I) 


I > LSN? 


CALL  COND(I) 
CALL  FCON(I) 
CALL  PVOL(I) 


KPR  = K 


CALL  TEST (I) 

\ 

I SAVE  = 07- 
yes 

I 

— < NS? 


NFIRES  = K 


CALL  ELEM 


CALL  CONDS(I) 
CALL  FCONS(I) 
CALL  PVOLS(I) 


KPR  = K 


m—m 


. 


I i 


Comment  No. 
23. 


24. 

25. 

26. 

27. 


28. 


29. 

30. 

31. 


32. 


33. 

34. 

35 


36. 


MAIN  PROGRAM'  FLOW  CHART  (Continued) 


I 


ISW  * 1 

1 


©- 


CALL  AFP 

\ 

•CALL  ATMOS 


I 


NPASS  « NPASS  + 1 

NPASS^  =*  1? E25- 

no 

I 


IPRl  = MOD ( ITIME , IPENS ) 
IPP.2  = MOD  (ITIME,  IPSPR) 


i 


CALL  OUTPUT 

\ 

CALL  PLOTS 

\ 


i ves  'NORMAL 

ITIME  > ITFIN?  * SPRINT: 


no 


COMPLETION ' 


no 


i 

IM 

(5) 


ITIME  ■ ITIME  + IDELT 

JL 

F 
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COMMENTS  FOR  PROGRAM  FLOW  CHART 


No. 

1.  ' NPASS ' is  a counter,  initialized  at  this  point,  which 
will  contain  the  number  of  times  the  cabin  atmosphere 
computations  have  been  performed.  IPR1,  IPR2  are 
print  controls. 

2.  Subroutine  'INPUTG'  initializes  and  defines  those 
variables  pertaining  to  the  geometry  of  the  cabin  section. 

3.  Test  for  run  termination  tno  additional  input  cards  to 
be  read) . 

4.  Eighty  characters  of  run  identification  are  printed. 

5.  Subroutine  ' INPUTM'  reads  all  input  data  pertaining  to 
the  material  properties  of  each  surface. 

6.  Subroutine  'INIT'  and  'INIT2'  perform  basic  computations 
from  the  input  data  and  initialize  those  variables 
required  to  start  the  integration. 

7.  Subroutine  ' INPUTO'  reads  all  data  relating  to  the 
ignition  source. 

8.  Subroutine  'PRTINP'  provides  an  input-echo  of  selected 
input  data. 

9.  This  point  is  the  start  of  the  primary  integration  loop 
in  the  program.  'ITIM2'  is  the  time  associated  with  the 
flame  propagation  computations,  ' ITIME f is  the  time 
associated  with  the  cabin  atmosphere  computations. 

10.  Initialize  sums  which  contain  emission  rates  for  all  fires. 

11.  'I*  is  the  surface  index  Cl  through  20  are  cabin  lining 
surfaces,  21  through  29  are  seat  groups). 

12.  Test:  Any  flaming  elements  on  surface  ’I’?  If  not,  by- 
pass flame  propagation  computations  for  this  surface. 

13.  Subroutine  ' FIRE ' isolates  a fire  on  the  specified 
surface  and  performs  the  computations  of  the  flame 
properties. 

14.  If  K * KPR  at  this  point,  a new  fire  was  not  found  in 
subroutine  'FIRE'. 

15.  Subroutine  'RATES'  determines  the  heat  flux  at  various 
points  associated  with  one  specific  fire  and  interpolates 
for  material  properties  as  a function  of  this  heat  flux. 

I 
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COMMENTS  FOR  PROGRAM  FLOW  CHART  (continued) 

No. 

16.  If  I LSN,  the  fire  under  consideration  is  located  on 
a seat;  otherwise  the  fire  is  located  on  a cabin  lining 
surface. 

17.  Subroutines  'COND',  'CONDS'  determine  flame  propagation 
via  conduction  (creeping  flame  spread) . 

Subroutine  'FCON',  'FCONS'  determine  flame  propagation 
via  flame  contact. 

Subroutines  'PV0L1,  'PVOLS'  test  for  possible  elemental 
change  of  state  due  to  the  pyrolysis  (smoldering)  of 
elements  in  the  vicinity  of  a fire. 

18.  Subroutine  'TEST*  determines  if  any  flaming  elements 
change  to  the  charred  state  and  sums  the  emission  rates 
for  each  fire. 

19.  If  ' ISAVE ' / 0;  then  return  control  to  subroutine  'FIRE' 
(continue  to  search  surface  ' I*  for  fires). 

20.  Test;  Have  all  cabin  lining  surfaces  and  seat  groups 
been  examined  during  this  time  step?  If  not,  increment 
for  next  surface. 

21.  The  variable  'NFIRES'  contains  the  total  number  of 
distinct  fires  in  progress  during  this  time  interval. 

22.  Subroutine  'ELEM  updates  the  time  counters  and  indicators 
associated  with  each  element. 

23.  ' ISW'  is  a switch.  When  ISW  = 1,  indicates  that  the 
flame  propagation  computations  have  been  performed 
this  time  step. 

24.  Subroutine  'AFP*  determines  the  total  number  of  flaming 
and  smoldering  elements  and  sums  emission  rates. 

25.  Subroutine  'ATMOS'  contains  all  of  the  equations  describing 
the  cabin  atmosphere. 

26.  Add  'one*  to  the  pass  counter. 

27.  Test:  If  this  is  the  first  pass  through  the  program, 
automatically  print  flame  propagation  and  cabin  atmosphere 
data. 

28.  Determine  if  flame  propagation  and/or  cabin  atmosphere 
data  is  to  be  printed  this  time  step. 
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COMMENTS  FOR  PROGRAM  FLOW  CHART  (continued) 


No. 

29.  Subroutine  'OUTPUT'  consists  of  the  required  print  and 
format  statements  and  controls  to  obtain  the  output 
data  as  required. 

30.  Subroutine  'PLOTS'  writes  selected  output  variables  to 
a file  for  later  plotting. 

31.  Test:  If  simulation  time  has  expired,  print  appropriate 
message  and  terminate  the  run. 

32.-34.  If  the  flame  propagation  computations  have  been  performed 
this  pass,  reset  computer  words  containing  the  element 
states  information. 


35.  If  flame  propagation  computations  are  to  be  performed 
the  next  time  step,  re-enter  appropriate  loop. 

36.  If  flame  propagation  computations  are  not  required  next 
time  step,  increment  cabin  atmosphere  'At'  and  re-enter 
cabin  atmosphere  computations. 
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C. 3 INPUT  DATA  PREPARATION 


This  section  describes  the  input  requirements  of  DACFIR2. 
The  preparation  of  each  input  card  is  described,  and,  where 
necessary,  a brief  explanation  of  the  input  data  requirements 
and  options  is  included.  Following  the  input  preparation 
instructions  is  a listing  of  a sample  input  data  deck.  Numbers 
in  the  left  column  of  the  data  description  identify  specific 
cards  in  the  sample  input  deck.  The  specific  set  of  input  data 
shown  in  the  listing  was  used  to  create  Case  26P. 

In  the  data  description  shown  below,  three  format  types 
are  referenced.  They  are  as  follows. 

Type  Description 

A Alphanumeric,  any  combination  of 

letters,  numbers,  and  special 
characters  (including  blanks) 
may  be  entered  in  the  appropriate 
column. 

I Integer,  the  entry  must  be  right 

justified  in  the  field  (range  of 
columns).  Example:  when  the 
number  '25'  is  entered  in  a five- 
column  field,  it  must  be  pre- 
ceded by  three  blanks. 

F Floating  point,  the  entry  may 

appear  anywhere  in  the  specified 
field,  but  the  insertion  of  a 
decimal  point  is  mandatory. 
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DACFIR2 

INPUT  DATA  DESCRIPTION 


Card 

Number 

Card 

Type 

1 

Variable 

Name 

Dimension 

Column 

Fmt 

Type 

Description 

1 

IDENT 

(20) 

1-80 

A 

Run  I . 0 . 

2 

2 

LSN 

— 

1-5 

I 

Number  lining  surfaces 

NSG 

— 

6-10 

I 

Number  seat  groups 

ICLL 

— 

11-15 

I 

Surface  number  of 
leftmost  ceiling  surface 

ICIR 

— 

16-20 

I 

Surface  number  of 
rightmost  ceiling  surface 

NTXG 

— 

21-25 

I 

Number  of  gases 
(excluding  0%)  to  be  used 

ID 

26-30 

I 

Flag  to  specify  the  use 
of  "square"  or  "round" 
cabin  cross  section  for 
gas-dynamics  calculations, 
ID-0  indicates  "round” 
case  to  be  used.  ID«1 
the  "souare." 

3 

3 

RFWS 

— 

11-20 

F 

Flame  spread  rate 
sidewall  to  seat  or 
seat  to  sidewall  (ft/ sec) 

DWS 

— — 

21-30 

F 

Separation  distance 
outboard  seats  to 
sidewall  (ft) 

CH 

— 

31-40 

F 

Cabin  Floor  to 
celling  height  (ft) 

CL 

— 

41-50 

F 

Cabin  section  length 
(see  Figure  2-1) 

CU 

— — 

51-60 

F 

Cabin  width  at  floor  (ft) 

SL 

— 

61-70 

F 

Oetailed  section 
length  (see  Figure  2-1 ) 

(Ft) 
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Card 

Number 

Card 

Type 

4 

Variable 

Name 

01 mens  ion 

Col umn  Fmt 

Type 

Description 

4-23 

SUO 

20 

1-10 

F 

Surface  "width"  - surface  dimension 
in  feet,  perpendicular  to  the 
y (down  cabin  direction)  - this  is 
the  actual  width  for  horizontal 
surfaces  and  the  height  for  vertical 
surfaces,  (ft) 

Z 

20 

11-20 

F 

z displacement  (height)  of  the 
lowest  part  of  surface  from  floor  (ft) 

VN 

20,3 

21-30 

31-40 

41-50 

F 

F 

F 

x component  of  surface  normal  (ft) 
y component  of  surface  normal  (ft) 
z component  of  surface  normal  (ft) 

IMATl 

20 

51-55 

I 

Surface  material  ident.  number 
‘Material  identification  number  is 
determined  by  the  order  of  input  of 
the  materials  data.  Thus  IMATL(i)*j 
indicates  that  surface  i is  of  material 
type  j and  this  material's  data  is  the 
jth  in  the  input  order  (see  below) 

24-29 

5 

SGWO 

9 

1-10 

F 

Seat  group  width  (ft) 

XCOR 

9 

11-20 

F 

x-coordinate  of  lefthand  forward 
comer  of  seat  group  (ft) 

rcoR 

9 

21-30 

F 

y-coordinate  of  lefthand  forward 
comer  of  seat  group  (ft) 

30 

6 

IMATS 

7 

i-2  | 

3-4 

13-14  ' 

I 

I 

i 

Seat  surfaces  materia! 
ident.  number  - follows 
format  of  lining  surfaces 
material  identification 

31 

7 

NV 

1-5 

i 

Number  of  openings  (vents)  in 
cabin  section  to  exterior  or 
other  cabin  volumes.  Includes 
doors  and  emergency  exits. 

32-33 

8 

VENTT 

10 

1-10 

F 

Distance  to  top  of  opening 
from  floor  ( ft) 

VENTH 

10 

11-20 

F 

Height  (z  dimension)  of 
opening  (ft) 

VENTVI 

10 

21-30 

F 

Width  (y  or  x dimension)  of 
opening  (ft) 

FLOW 

10 

31-40 

F 

Imposed  air  flow  rate  through 
opening;  may  be  positive,  negative, 
or  zero  (ft’ /min) 

013 
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Card 

Card 

Variable 

01 mens ion 

Column 

Fmt 

Description 

Number 

TyP* 

9 

Name 

Type 

34-40 

QTAB 

7 

i-10 

F 

Heat  of  combustion  (BTU/lbm) 
of  material 

GTAB 

7 

11-20 

F 

Stoichiometric  fuel  to  oxygen 
ratio  of  material  (no  units) 

RTAB 

7 

21-30 

F 

Fuel  vapor  (pyrolyzate)  density 
at  burning  material  surface 
(lbm/ft3) 

UTAB 

7 

31-40 

F 

Fuel  vapor  (pyrolyzate) 
flow  velocity  at  burning 
material  surface  (ft/sec) 

RAOTAB 

7 

41-50 

F 

Fraction  of  material's  heat  of 
combustion  released  to  envirms 

* One 

card  of  type  9 must  be  prepared 

for  each  material. 

by  flame  radiation  (no  units) 

41 

10 

TP 

7 

1-10 

F 

Time  of  transition  of  an 

11-20 

F 

element  from  the  original 

I 

state  to  the  smoldering  state 

61-70 

F 

(sec).  Enter  values  in  the 
order  of  material  identification 

42 

11 

TPC 

7 

1-10 

F 

Time  of  transition  of  an  element 

11-20 

F 

from  the  smoldering  to  the 

: 

* 

charred  state  (sec).  Enter 

61-70 

F 

values  in  the  order  of  materials 
identification 

43 

12 

RSS 

7 

1-10 

F 

Smoke  production  rate  for 

.* 

11-20 

F 

each  material  In  the 

! 

• 

smoldering  state 

61-70 

F 

(part1cles/ft2-sec)  Enter 
values  in  the  order  of  materials 
i .d. 

44-50 

13 

RSG 

10,7 

1-10 

F 

Toxic  gas  production  rate 

11-20 

F 

for  each  material  in  the 

61-70 

F 

smoldering  state 
(10*6  ibm/ft2-sec) 

One  card  of  this  type  for  each  toxic  gas. 
Rate  data  for  all  materials  for  the  1th 

gas  on 

a single  card  in  the  order  of 

materials  Identification. 

51-57 

14 

TABX 

18,7,6 

1-5 

F 

Radiant  flux  (BTU/ft2-sec) 

TABY 

18,7,6 

6-13 

F 

i Material  horizontal  flame  spread 
| rate  at  above  flux  level  (ft/sec) 

TABX 

18,7,6 

14-18 

F 

Radiant  flux  (BTU/ft2-sec) 

TABY 

18,7,6 

19-26 

F 

| Flame  spread  rate  at  above  flux 

I (ft/ sec) 

This  format 

followed 

for  six  pairs 

per  card 
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Card  Card  Variable  Dimension  Column  Fmt  Description 

Number  Type  Name  Type  


58-64 

15  - format  identical  to  card  type  14  - 

Vertical  upward  flame  spread  rate  (ft/sec) 

65-71 

16  - format  identical  to  card  type  14  - 

Vertical  downward  flame  spread  rate  (ft/sec) 

72-78 

17  - format  identical  to  card  type  14  - 

Time  to  iqnite  (virqin  to  flaming)  (sec) 

79-85 

18  - format  identical  to  card  type  14  - 

Heat  release  rate-flaminq  state  (8TD/ft2-sec) 

36-92 

19  - format  identical  to  card  type  14  - 

Smoke  release  rate-flaming  (Part1cles/ft2-sec) 

93-99 

20  - format  identical  to  card  type  14  - 

Time-pyrolysis  to  extinction  (sec) 

(No  data  for  this  item,  zero  entered) 

100-106 

21  - format  identical  to  card  type  14  - 

Time  to  bum  out  (flaming  to  charred)  (sec) 

107-113 

22  - format  identical  to  card  type  14  - 

Release  rate  (flamlnq)  for  1st  toxic  gas  (10'®  lbm/ft2-sec) 

114-120 

23  - format  Identical  to  card  type  14  - 

Release  rate  (flaminq)  for  2nd  toxic  gas  (10*®  ibm/ft2-sec) 

121-127 

24  - format  Identical  to  card  type  14  - 

Release  rata  (flamlnq)  for  3rd  toxic  qas  (10'®  lbm/ft‘-sec) 

128-34 

25  - format  identical  to  card  type  14  - 

Release  rate  (flamlnq)  for  4th  toxic  gas  (10*6  lbm/ftZ-sec) 

133-141 

• 

26  - format  identical  to  card  type  14  - 2 . 

Release  rata  (flaming)  for  5th  toxic  das  (10-6  lbm/ft  -sec) 

142-148 

27  - format  Identical  to  card  type  14-  - , 

Release  rate  (flaming)  for  6th  toxic  gas  (10  lbm/ft  -sec) 

149-155 

28  - format  Identical  to  card  type  14-  - - 

Release  rate  (flaming)  for  7th  toxic  gas  (10*°  lbm/ft  -sec) 

•NOTE: 

The  sample  deck  contains  data  on  7 toxic  gases:  the  procedure 
above  can  be  continued  for  up  to  10  gases. 
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Card 

Number 

Card 

Type 

* 

32 

i! 

<0 

> 

Dimension 

Column  Fmt 

Type 

Description 

156-162 

NGAS 

10 

1-6  A 

Alphanumeric  name  of  toxic 
gas,  max.  6 characters. 

XMW 

11-22 

E 

Gas  molecular  weight  (Ibm/lbraole) 

GIVI 

23-24 

E 

First  check  level  (PPM) 

GLV2 

35-46 

E 

Second  check  level  (PPM) 

GLV3 

47-58 

E 

Third  check  level  (PPM) 

•NOTE: 

Up  to  10  cards  of  type  32  can 
must  correspond  to  the  order 
has  been  entered  previously. 

be  entered.  The  order  In  which  they  are  entered 

In  which  the  release  data  (flaming  and  smoldering) 

163 

33 

OP 

7 

1-10  F 

11-20  F 

61-70  F 

Threshold  flux  levels 
for  transition  to 
smoldering  state  (BTU/ft2-sec) 
for  each  material 

164 

34 

CPU 

— 

1-10  F 

Specific  heat  of  materials  at 
ambient  conditions,  average 
value,  BTU/( lbm  ■ °R) 

PHOM 

— 

21-30  F 

Bulk  density  of  materials, 
average  value(lbm/ft3) 

XK 

-— - 

21-30  F 

Thermal  conductivity  of 
materials  at  ambient  conditions, 
average  value,  8TU/(ft  • sec  • OR) 

XPEN 

— 

31-40  F 

Heat  penetration  depth  of 
materials,  average  value,  ft 

TO 

— 

41-50  F 

Ambient  temperature,  OR 

165 

35 

OELTAT 

TFINAL 

— 

1-10  F 

11-20  F 

Integration  step  size  (small  step) 
for  gas  dynamics  model  (sec) 
Minimum  step  size  is  0.001  sec. 
Total  time  for  simulation  (sec) 
(starting  time  assumed  to  be 

0.  sec  always) 

I RATI 0 

11-15  I 

Ratio  of  steps  through  gas 
dynamics  calculation  to  steps 
through  flame  spread  calculation  - 
must  always  be  greater  than  or 
egual  to  one 

IPEMS 

— — 

15-20  1 

Output  printing  Interval  for 
cabin  atmosphere  sunnary  expressed 
as  a multiple  of  DELTAT 

IPSPR 

21-25  I 

Output  printing  Interval  for  flame 
spread-fire  growth  summary 
expressed  as  a multiple  of  DELTAT 
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Card 

Number 

Card 

Type 

Variable 

Name 

Dimension  Column 

Fmt 

IYJ22. 

Description 

166 

36 

Q8KGN0 

1-10 

F 

Externally  imposed  "background" 
radiant  flux  level  (8TU/ft2-sec) 

167 

37 

QCI 

1-10 

F 

Heat  of  combustion  of  ignition 
source  fuel  (BTU/lbm) 

SAMI 

11-20 

F 

Stoichiometric  fuel  to  oxygen 
ratio  for  ignition  source  fuel 
(no  units) 

mi 

21-30 

F 

Ignition  fuel  vapor  density  at 
burning  fuel  surface  (Lbm/ft'3) 

XMUI 

31-40 

F 

Ignition  fuel  vapor  velocity 
at  burning  fuel  surface  (ft/sec) 

RAO  I 

41-50 

F 

Fraction  of  ignition  fuel's  heat 
of  combustion  released  to  environs 
by  flame  radiation  (no  units) 

XMFI 

51-60 

F 

•Total  amount  of  ignition  source 
fuel  (lbm) 

•Only  one  ignition  fire  Is  allowed; 
fully  developed  at  time  ■ 0. 

the  ignition  fire  is  assumed  to  be 

168 

. 

38 

001 

1-10 

F 

Heat  release  rate  for  ignition 
source  fuel  (BTU/ft2-sec) 

169 

39 

RSI 

1-10 

F 

Smoke  release  rate  for  ignition 
source  fuel  (particles/ft2-sec) 

Mot  shown 

39 

RTGI 

RTGI 

10  11-20 

21-30 
31-40 

71-80 

10  1-10 

11-20 

F 

F 

F 

F 

F 

Toxic  gas  release  rates 
(10*6  lbei/ft2-sec)  for 
ignition  source  fuel.  If 
there  are  more  than  eignt 
toxic  gases  the  last  one  or 
two  are  entered  in  the 
first  two  fields  of  a second 
card 

170 

40 

I GSM 

1-5 

I 

Surface  number  uoon  which  the 
ignition  source  fuel  is  located 

171 

41 

NIJSQ 

— — 1-5 

I 

Number  of  square  elements 
(0.5  ft  x 0.5  ft)  which  make  up 
the  base  of  the  ignition  source 
fire 

PIGN 

6-15 

F 

The  length  of  the  perimeter  of  the 
base  of  the  ignition  source  fire, 
(ft) 
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Card 

Number 

Card  Variable 

Type  Name 

Dimension  Column  Fmt 

Type 

Description 

172-177 

42  1 

1-5 

I 

1 Index  of  Ignition  source 
fire  base  element 

J 

6-10 

I 

j Index  of  ignition  source 
fire  base  element 

Enter  one  pair  of  (i.J)  Indices  per 
will  be  equal  to  the  value  of  NIGSQ 
Figure  . 

card.  The  total  number  of  cards 

For  number  of  elements  see 

178 

43  NIJC 

1-5 

I 

Total  number  of  elements,  on 
any  or  all  surfaces,  which  are 
to  be  set  to  the  charred  (Inert) 
state  at  the  start  of  the 
simulation.  If  the  value  is 
entered  as  zero  no  cards  of  type 

44  are  included  in  the  deck 

Not  shown 

a I 

1-5 

I 

1 Index  of  charred  (Inert)  element 

J 

6-10 

I 

j index  of  charred  (inert)  element 

Enter  one  pair  of  (i,j)  Indices  per 
will  be  equal  to  the  value  of  NIJC . 
Figure 

card.  The  total  number  of  cards 

For  numbering  of  elements  see 

DACFIR2  - Sample  Input  Data  Deck 


Card  Column 


NO. 

l 

11 

21 

31 

41 

si 

61 

71 

I 

26  FT  CABIN  MOCK- 

UP  FIRE 

TEST  - IN-SERVICE  (1968) 

MATERIALS  (AIA-C3P 

co 

v£> 

1 

GO 

• 

ro 

2 

20 

6 13 

9 7 

0 

3 

0.2 

0.1 

7.5 

25.83 

11.5 

7.5 

4 

11.5 

0. 

0. 

0. 

1.0 

1 

5 

1.5 

0. 

1 .0 

0. 

0. 

2 

6 

3.0 

0. 

-1.0 

0. 

0. 

2 

7 

1.5 

0. 

-1.0 

0. 

0. 

2 

8 

2.0 

6.0 

0. 

0. 

-1.0 

3 

9 

0.5 

0. 

-1.0 

0. 

0. 

4 

10 

2.0 

6.5 

0. 

0. 

1.0 

5 

11 

0.5 

0. 

-1.0 

0. 

0. 

6 

12 

1.0 

7.0 

0. 

0. 

-1.0 

6 

13 

0.5 

0. 

-1.0 

0. 

0. 

6 

14 

9.5 

7.5 

0. 

0. 

-1.0 

6 

15 

0.5 

0. 

1.0 

0. 

0. 

6 

18 

1.0 

7.0 

0. 

0. 

-1.0 

6 

17 

0.5 

0. 

1.0 

0. 

0. 

6 

18 

2.0 

6.5 

0. 

0. 

1.0 

5 

19 

0.5 

0. 

1.00 

0. 

0. 

4 

20 

2.0 

6.0 

0. 

0. 

-1.0 

3 

21 

1.5 

0.  . 

1.00 

0. 

0. 

2 

22 

3.0 

0. 

1.0 

0. 

0. 

2 

23 

1.5 

0. 

1.00 

0. 

0. 

2 

24 

3.5 

0. 

0. 

25 

5.0 

6.5 

0. 

26 

3.5 

0. 

2.5 

27 

5.0 

6.5 

2.5 

28 

3.5 

0. 

5.0 

29 

5.0 

6.5 

5.0 

30 

7 7 7 7 

7 7 7 

31 

2 

32 

7.0 

7.0 

3.0 

0.0 

33 

4.0 

4.0 

2.0 

0.0 

34 

7000. 

2. 

0.035 

0.25 

0.25 

35 

7000. 

2. 

0.035 

0.25 

0.25 

36 

7000. 

2. 

0.035 

0.25 

0.25 

37 

7000. 

2. 

0.035 

0.25 

0.25 

38 

7000. 

2. 

0.035 

0.25 

0.25 

39 

7000. 

2. 

0.035 

0.25 

0.25 

40 

7000. 

2. 

0.035 

0.25 

0.25 

41 

59.0 

49.0 

75.2 

21.0 

3.0 

41.0 

21.0 

42 

455.0 

163.0 

662.0 

334.0 

113.0 

136.0 

162. 

43 

2.57 

5.8 

74.6 

17.1 

23.1 

3.23 

10.8 

44 

1457. 

0.0 

0.0 

37.2 

72.1 

256.6 

673.0 

45 

41.2 

5.1 

.419 

6.23 

75.0 

61.2 

145.6 

46 

27.8 

88.8 

0.0 

27.7 

58.4 

39.9 

0.035 

47 

1.18 

0.0 

0.0 

.80 

.432 

.74 

1.65 

48 

0.0 

8.5 

0.0 

0.0 

0.0 

0.0 

0.0 

49 

.55 

0.0 

0.0 

0.0 

0.0 

0.0 

0.63 

50 

0.0 

0.0 

1.92 

0.0 

0.0 

0.0 

0.0 

51 

0.0  0.005  1.32 

0.0082 

2.20  0.0079 

3.08  0.0198 

4.4 

0.0144  6.0 

0.0144 

52 

0.0  0.0 

1.32 

0.01 

2.20  .0106 

3.08  .0408 

4.4 

.0632  6.0 

.0632 

53 

0.0  0.0 

1.32 

0.0 

2.20  0.0033 

3.08  0.0079 

4.4 

0.0145  6.0 

0.0145 

54 

0.0  .001 

1.32 

.0038 

2.20  .0076 

3.08  .0326 

4.4 

.0216  6.0 

.0216 

55 

0.0  0.1 

1.32 

.1042 

2.20  .1389 

3.08  .1389 

3.96 

.4167  6.0 

.4167 
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Card  Column 


NO. 

i 

u 

21 

31 

41 

si 

61 

71 

56 

0.0 

0.005 

1.32 

0.0139 

2.20 

.0463 

3.08 

.0556 

3.96 

.0740 

6.0 

.0740 

57 

0.0 

0.0179 

1.32 

0.0179 

2.2 

0.0246 

3.08 

0.0463 

4.4 

0.0716 

6.0 

0.0716 

58 

0.0 

0.005 

1.32 

0.0082 

2.20 

0.0079 

3.08 

0.0198 

4.4 

0.0144 

6.0 

0.0144 

59 

0.0 

0.0 

1.32 

.0111 

2.20 

.0833 

3.08 

.0764 

4.4 

.0875 

6.0 

.0875 

60 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0098 

3.08 

0.0201 

4.4 

0.0471 

6.0 

0.0471 

61 

0.0 

.01 

1.32 

.0227 

2.20 

.0179 

3.08 

.0336 

4.4 

.0398 

6.0 

.0398 

62 

0.0 

0.1 

1.32 

.1042 

2.20 

.1389 

3.08 

.1389 

3.96 

.4167 

6.0 

.4167 

63 

0.0 

0.005 

1.32 

0.0139 

2.20 

.0463 

3.08 

.0556 

3.96 

.0740 

6.0 

.0740 

64 

0.0 

0.0179 

1.32 

0.0179 

2.2 

0.0417 

3.08 

0.105 

4.4 

0.333 

6.0 

0.333 

65 

0.0 

0.005 

1.32 

0.0082 

2.20 

0.0079 

3.08 

0.0198 

4.4 

0.0144 

6.0 

0.0144 

66 

0.0 

0.0 

1.32 

.0038 

2.20 

.0227 

3.08 

.0321 

4.4 

.0820 

6.0 

.0820 

67 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0015 

3.08 

0.0087 

4.4 

0.0146 

6.0 

0.0146 

68 

0.0 

.002 

1.32 

.0025 

2.20 

.0064 

3.08 

.01 

4.4 

.0419 

6.0 

.0419 

69 

0.0 

0.1 

1.32 

.1042 

2.20 

.1389 

3.08 

.1389 

3.96 

.4167 

6.0 

.4167 

70 

0.0 

0.005 

1.32 

0.0139 

2.20 

.0463 

3.08 

.0556 

3.96 

.0740 

6.0 

.0740 

71 

0.0 

0.00633 

1.32 

0.00633 

2.2 

0.0330 

3.08 

0.0590 

4.4 

0.264 

6.0 

0.264 

72 

0.0 

20.0 

1.32 

12.0 

2.20 

2.0 

3.08 

1.3 

4.4 

2.30 

6.0 

2.30 

73 

0.0 

30.0 

1.32 

26.7 

2.20 

20.00 

3.08 

11.56 

4.4 

9.89 

6.0 

9.89 

74 

0.0 

30.0 

1.32 

20.0 

2.20 

16.25 

3.08 

14.75 

4.4 

35.33 

6.0 

35.33 

75 

0.0 

20.0 

1.32 

14.3 

2.20 

5.0 

3.08 

8.7 

4.4 

3.9 

6.0 

3.9 

76 

0.0 

15.0 

1.32 

15.0 

2.20 

1.0 

3.08 

4.0 

3.96 

1.0 

6.0 

1.0 

77 

0.0 

8.0 

1.32 

6.0 

2.20 

2.0 

3.08 

4.0 

3.96 

1.0 

6.0 

1.0 

78 

0.0 

20.0 

1.32 

16.3 

2.2 

5.00 

3.08 

7.19 

4.4 

4.58 

6.0 

4.58 

79 

0.0 

3.0 

1.32 

4.0 

2.20 

5.18 

3.08 

6.28 

3.96 

0.9387 

6.0 

0.9387 

80 

0.0 

.20 

1.32 

.433 

2.20 

2.37 

3.08 

3.15 

4.4 

2.73 

6.0 

2.73 

81 

0.0 

0.0 

1.32 

0.0 

2.20 

0.917 

3.08 

2.84 

4.4 

3.51 

6.0 

3.51 

82 

0.0 

5.0 

1.32 

7.47 

2.20 

9.77 

3.08 

13.6 

4.4 

7.3 

6.0 

7.3 

83 

0.0 

2.0 

1.32 

2.871 

2.20 

4.051 

3.08 

2.896 

3.96 

3.136 

6.0 

3.136 

84 

0.0 

0.5 

1.32 

1.694 

2.20 

2.389 

3.08 

3.167 

3.96 

3.234 

6.0 

3.234 

85 

0.0 

5.0 

1.32 

5.37 

2.2 

5.50 

3.08 

5.46 

4.4 

4.47 

6.0 

4.47 

86 

0.0 

10.0 

1.32 

20.0 

2.20 

26.84 

3.08 

38.92 

3.96 

83.05 

6.0 

83.05 

87 

0.0 

2.0 

1.32 

5.92 

2.20 

15.28 

3.08 

37.95 

4.4 

67.0 

6.0 

67.0 

88 

0.0 

0.0 

1.32 

0.0 

2.20 

2.41 

3.08 

33.9 

4.4 

55.5 

6.0 

55.5 

39 

0.0 

98.0 

1.32 

98.4 

2.20 

98.2 

3.08 

426.0 

4.4 

116.2 

6.0 

116.2 

90 

0.0 

5.0 

1.32 

11.66 

2.20 

0.0 

3.08 

14.07 

3.96 

19.46 

6.0 

19.46 

91 

0.0 

10.0 

1.32 

8.695 

2.20 

11.25 

3.08 

2.186 

3.96 

16.45 

6.0 

16.45 

92 

0.0 

4.0 

1.32 

5.68 

2.2 

8.02 

3.08 

10.9 

4.4 

20.6 

6.0 

20.6 

93 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

94 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

95 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

96 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

97 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

3.96 

0.0 

6.0 

0.0 

98 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

3.96 

0.0 

6.0 

0.0 

99 

0.0 

0.0 

1.32 

0.0 

2.2 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

100 

0.0 

750.0 

1.32 

615.0 

2.20 

327.0 

3.08 

279.0 

4.4 

186.3 

6.0 

136.3 

101 

0.0 

130.0 

1.32 

126.3 

2.20 

83.0 

3.08 

56.3 

4.4 

52.6 

6.0 

52.6 

102 

0.0 

20000 

1.32 

2000 

2.20 

1106 

3.08 

606.0 

4.4 

538.1 

6.0 

538.1 

103 

0.0 

250.0 

1.32 

221.3 

2.20 

211.3 

3.08 

152.0 

4.4 

134.6 

6.0 

134.6 

104 

0.0 

150 

1.32 

180.0 

2.20 

124.0 

3.08 

321.0 

3.96 

270.0 

6.0 

270.0 

105 

0.0 

150 

1.32 

174.0 

2.20 

245.0 

3.08 

225.0 

3.96 

324.7 

6.0 

324.7 

106 

0.0 

562. 

1.32 

551. 

2.2 

542. 

3.08 

607. 

4.4 

813. 

6.0 

813. 

107 

0.0 

500.0 

1.32 

680.8 

2.20 

2241.0 

3.08 

2071.0 

4.4 

4249.0 

6.0 

4249.0 

108 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

348.6 

4.4 

0.0 

6.0 

0.0 

109 

0.0 

0.0 

1.32 

0.0 

2.20 

25.6 

3.08 

0.0 

4.4 

178.0 

6.0 

178.0 

no 

0.0 

1500. 

1.32 

1500. 

2.20 

1548. 

3.08 

1535. 

4.4 

36.78 

6.0 

36.78 
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Card  Column 


No. 

l 

11 

21 

31 

<*i 

51 

61 

71 

Ill 

0.0 

150 

1.32 

200 

2.20 

231.6 

3.08 

636.6 

3.96 

545.7 

6.0 

545.7 

112 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

3.96 

0.0 

6.0 

0.0 

113 

0.0 

2500 

1.32 

2500. 

2.2 

3049. 

3.08 

1979. 

4.4 

2280. 

6.0 

2 280. 

114 

0.0 

15.0 

1.32 

19.63 

2.20 

56.02 

3.08 

87.54 

4.4 

170.6 

6.0 

170.6 

115 

0.0 

50.0 

1.32 

50.0 

2.20 

80.25 

3.08 

135.6 

4.4 

213.7 

6.0 

213.7 

116 

0.0 

0.0 

1.32 

0.0 

2.20 

5.02 

3.08 

15.2 

4.4 

53.9 

6.0 

53.9 

117 

0.0 

150. 

1.32 

150. 

2.20 

197.0 

3.08 

230.1 

4.4 

509.3 

6.0 

509.3 

118 

0.0 

100 

1.32 

150. 

2.20 

184.2 

3.08 

137.0 

3.96 

88.96 

6.0 

88.96 

119 

0.0 

20 

1.32 

36.69 

2.20 

121.8 

3.08 

141.2 

3.96 

113.0 

6.0 

113.0 

120 

0.0 

169. 

1.32 

169. 

2.2 

169. 

3.08 

182. 

4.4 

169. 

6.0 

169. 

121 

0.0 

7.5 

1.32 

11.77 

2.20 

33.19 

3.08 

45.39 

4.4 

77.68 

6.0 

77.68 

122 

0.0 

190.1 

0 1.32 

194.0 

2.20 

174.4 

3.08 

64.26 

4.4 

481.5 

6.0 

481.5 

123 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.060 

4.4 

0.202 

6.0 

0.202 

124 

0.0 

80. 

1.32 

80. 

2.20 

85.6 

3.08 

119.0 

4.4 

134.4 

6.0 

134.4 

125 

0.0 

50 

1.32 

50. 

2.20 

80.04 

3.08 

67.42 

3.96 

37.21 

6.0 

37.21 

126 

0.0 

25 

1.32 

51.98 

2.20 

59.07 

3.08 

80.40 

3.96 

90.0 

6.0 

90.0 

127 

0.0 

0. 

1.32 

0. 

2.2 

0. 

3.08 

0. 

4.4 

0.250 

6.0 

0.250 

128 

0.0 

G.O 

1.32 

0.2176 

2.20 

1.228 

3.08 

1.919 

4.4 

4.310 

6.0 

4.310 

129 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

130 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

131 

0.0 

0.0 

1.32 

5.0 

2.20 

13.9 

3.08 

6.2 

4.4 

4.4 

6.0 

4.4 

132 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

1.662 

3.96 

0.9175 

6.0 

0.9175 

133 

0.0 

0.0 

1.32 

0.0 

2.20 

.8739 

3.08 

.8326 

3.96 

.8490 

6.0 

.8490 

134 

0.0 

3.5 

1.32 

3.50 

2.2 

3.70 

3.08 

3.30 

4.4 

1.36 

6.0 

1.36 

135 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

136 

0.0 

30.0 

1.32 

30.0 

2-20 

35.83 

3.08 

52.82 

4.4 

37.69 

6.0 

37.69 

137 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

138 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

139 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

3.96 

0.0 

6.0 

0.0 

140 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

3.96 

0.0 

6.0 

0.0 

141 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

142 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

1.224 

6.0 

1.224 

143 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

144 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

145 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

146 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

3.96 

0.0 

6.0 

0.0 

147 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

3.96 

0.0 

6.0 

0.0 

148 

0.0 

0.0 

1.32 

0.0 

2.2 

0.0 

3.08 

0.0 

4.4 

0.631 

6.0 

0.631 

149 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

150 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

151 

0.0 

0.0 

1.32 

0.0 

2.20 

11.9 

3.08 

16.8 

4.4 

35.4 

6.0 

35.4 

152 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

153 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

3.96 

0.0 

6.0 

0.0 

154 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

3.96 

0.0 

6.0 

0.0 

155 

0.0 

0.0 

1.32 

0.0 

2.20 

0.0 

3.08 

0.0 

4.4 

0.0 

6.0 

0.0 

156 

CO  2 

44.0 

1.0E+3 

1.0E+4 

1.0E+5 

157 

CO 

28.0 

5.0E+1 

5.0E+1 

8.0E+3 

158 

HCL 

36.5 

5.0E+O 

3.5E+1 

1.0E+3 

159 

HCN 

27.0 

1.0E+1 

1.0E+1 

2.8E+2 

160 

HF 

20.0 

3.0E+0 

3.2E+1 

5.0E+1 

161 

SC  2 

46.0 

5.0E+0 

5.0E+1 

2.5E+2 

162 

SO  2 

48.0 

5.0E+0 

2.0E+1 

1.0E+2 

163 

2.2 

2.2 

2.2 

2 

.2 

2.2 

2 

.2 

2.2 

2.2 

164 

0.25 

35.0 

0.000084  0 

1.00833 

530. 

,0 

165 

0.005 

250. 

1000 

5 

5 
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Column 


Card 


No. 

1 

u 

166 

0. 

167 

18000. 

3.5 

168 

214.6 

169 

86.7 

18720 

170 

1 

171 

6 

5.0 

172 

22 

6 

173 

23 

6 

174 

22 

7 

175 

23 

7 

176 

22 

8 

177 

23 

8 

178 

0 

21 

31 

41 

0.271 

0.0198 

0. 

,40 

34.1 

0.0 

0. 

,0 

51  St  71 

10.8 

0.0  0.0  0.0 
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C.4  PROGRAM  OUTPUT 

The  following  pages  contain  a sample  of  the  printed  output 
of  DACFIR2 . This  output  was  produced  by  the  sample  input 
deck  shown  above.  On  pages  C-24,  C-25,  and  C-26  the  automatic 
input  summary  is  given.  All  input  data  except  the  tables  of 
materials'  properties  is  printed.  A plan  of  the  seat  location 
is  shown,  and  the  burning  time  of  the  ignition  fire,  computed 
from  the  input  data,  is  given.  Pages  C-27,  C-28,  and  C-29  show  the 
output  at  60  seconds  of  Case  26P. 

DACFIR2  also  writes  the  values  of  certain  variables  to 
a disk  file  or  tape  (as  determined  by  the  user)  for  off-line 
plotting  or  other  purposes.  The  variables  are  written  to 
unit  8 using  a formatted  FORTRAN  write  statement  with  an 
E12.5  format  for  each  variable.  The  variables  are  as  follows. 


Variable  and  Unit 


Record  1 

Time  (seconds) 

Upper  zone  gas 
temperature  (°F) 
Lower  zone  gas 
temperature  (°F) 
Smoke  concentration 
(optical  density/ft) 
Smoke  concentration 
(%  transmission/ft) 


Variable  Name 


TIME 

TUF 

TLF 

SOD 

PTRANS 


Record  2 


Time  (seconds)  TIME 

Oxygen  concentration  OPRCNT 

(%  vol.) 

Depth  of  the  lower  gas  HEIGHT 

zone  (ft) 

Upper  zone  materials  TSUF 

surface  temperature  (°F) 

Lower  zone  materials  TSLF 

surface  temperature  (°F) 


C-23 


LOOKING  FROM  THE  FRONT  OF  THE  CABIN  TOWARD  THE  REAR.THE  FLOOR  IS  SURFACE  NO  1, 
THE  LWR  RGT  SIDEWALL  IS  SURFACE  NO  2, THE  RGT  SIDE  CEILING  SURFACE  IS  9,  THE 
LEFT  SIDE  CEILING  SURFACE  IS  13,AND  THE  LWR  LEFT  SIDEWALL  IS  SURFACE  NO  20 

THE  EL£MDTCS(IN  THE  I DIRECTION) ARE  NUMBERED  Q)  THE  SANE  MAtHER  AS  THE  SURfACES 
THE  FRONT  MOST  ELEMENT  IS  J-l.THE  REAR  MOST  ELQSNT  IS  J-15 


SURFACE  DAZA(EXCL  SEATS) SURF  1 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

U 

12 

13 

14 

15 

16 

17 

18 

19 

20 


WIDTH, FT  UGT,FT  UNIT  NORM 


11. s 

0.0 

0.0 

0 

1.0 

1.5 

8.0 

-1.0 

0 

0.0 

3.0 

0.0 

-1.0 

0 

0 

0.0 

1.5 

0.0 

-1.0 

0 

0.0 

2.0 

6.0 

0.0 

0 

-1.0 

.5 

0.0 

-1.0 

0 

0.0 

2.0 

6.5 

0.0 

0 

1.0 

.5 

0.0 

-1.0 

0 

0.0 

1.0 

7.0 

0.0 

0 

-1.0 

.5 

0.0 

-1.0 

0 

0.0 

9.5 

7.5 

0.0 

0 

-1.0 

.5 

0.0 

1.0 

0 

0.0 

1.0 

7.0 

0.0 

0 

-1.0 

.5 

0.0 

1.0 

0 

0.0 

2.0 

6.5 

0.0 

0 

1.0 

.5 

0.0 

1.0 

0 

0.0 

2.0 

6.0 

0.0 

0 

-1.0 

1.5 

1.0 

1.0 

0 

0.0 

3.0 

0.0 

1.0 

0 

0.0 

1.5 

0.0 

1.0 

0 

0.0 

MAIL  TWE 

ixm 

XMAX 

xw 

XMX 

1 

l 

23 

0.0 

0.0 

2 

24 

26 

0.0 

1.5 

2 

27 

32 

1.5 

4.5 

2 

33 

35 

4.5 

6.0 

3 

36 

39 

6.0 

6.0 

4 

40 

40 

6.0 

6.5 

5 

41 

44 

6.5 

6.5 

6 

45 

45 

6.5 

7.0 

6 

46 

47 

7.0 

7.0 

6 

48 

48 

7.0 

7.5 

6 

49 

67 

7.5 

7.5 

6 

68 

68 

7.0 

7.5 

6 

69 

70 

7.0 

7.0 

6 

71 

71 

6.5 

7.0 

5 

72 

75 

6.5 

6.5 

4 

76 

76 

6.0 

6.5 

3 

77 

88 

6.0 

6.0 

2 

81 

83 

4.5 

6.0 

2 

84 

89 

1.5 

4.5 

2 

9# 

92 

0.0 

1.5 

THIS  PACE  IS  BEST  QUALITY  PRACTICABLE 
FROM  COPY  FURNISHED  TO  DDC  


SEAT  GROUP  DATA— 


ST  GAP  NO 

1 

2 

3 

4 

5 

6 


WIDTH, FT 
3.5 
5.0 
3.5 
5.0 
3.5 
5.0 


LOCATION 

XOOR.FT  TOOR.FT 


0.0 

6.5 

0.0 

6.5 

0.0 

6.5 


0.0 

0.0 

2.5 

2.5 

5.0 

5.0 


EACH  SEAT  CROUP  HAS  7 SURFACES: 


SURF  NO  SURFACE  MAIL  TSPE 

1 CUSHION  BOTTOM  7 

2 US  REAR  BACKREST  7 

3 OPR  REAR  BACKREST  7 

4 TOP  OF  BACKREST  7 

5 FRONT  OF  BACKREST  7 

6 CUSHION  TOP  7 

7 CUSHION  FRONT  7 

OVERHEAD  VIEW  OF  SEATS  SUPERIMPOSE)  ON  FLOOR  (NON-ZERO  INTEGER  IS  SEAT  GROUP  NUMBER) 

55555550000006666666666 
55555550000006666666666 


14 

13 

12 

U 

10 

9 

3 
7 
6 
5 

4 
3 
2 
1 


55555550000006666666666 

33333330000004444444444 

33333330000004444444444 

33333330000004444444444 

33333330000004444444444 

00000000000000000000000 

11111110000002222222222 

11111110000002222222222 

11111110000002222222222 

11111110060002222222222 

FRONT 


NATOS  DATA — 


NATO,  NO 

HEATS  VALUE 

STOCH  RATIO 

PYRQLZAIE  DENS 

PYROLZATE  VEL 

(BTU/LB) 

(LaM/CU  FT, 

(FT/SEC) 

1 

7000.0 

2.00 

.0350 

.25 

2 

7000.0 

2.00 

.0350 

.25 

3 

7000.0 

2.00 

.0350 

.25 

4 

7000.0 

2.00 

.0350 

.25 

5 

7000.0 

2.00 

.0350 

.25 

6 

7000.0 

2.00 

.0350 

.25 

7 

7000.0 

2.00 

.0350 

.25 

SPEC  HEAT (BTU/LBM-R) - .250 
MBIT*  (LBM/CH  FT) -35.00 
THERM  COND (BTO/TT-SEC-R) ■ .000084 


THERFAL  PENETRATION  DEPTH  (FT)-  .00833 

Aieiarr  tempcdegr  ri  -530.0 


jbispage 

from  COF^ 


TIME"  60.000  SEC  AFTER  IGNITICN 


CABIN  ADOSPHERE  SWWARY 


STRATI  FIS 

ZCNE  DEPTH  CAS  DENSITY 

CAS  TEMP  MAIL  SURF  TEMP 

HEAT  RATE  10  SURF 

CAS  .MODEL 

(FT)  (LBM/CD  FT) 

(DEC  F) 

(DEC  F) 

(BTU/SQ  FT-SEC) 

(JPPER  ZONE 

5.617  .0442 

438.6 

90.9 

.425 

LONER  TONE 

1.883  .0748 

70.9 

79.8 

.206 

TOTAL  FLOW  RATE  IBRD  VENTS 

OBC  CAS 

FLOW, ALL  FIRES 

(LBM/SEC) 

(LBM/SEC) 

OUT 

.373 

.804 

IN  0.000 

SMOKE  CONCENTRATION  IN  UPPER  ZCNE  TOXIC  CAS  CONCENTRATION  IN  UPPER  ZCNE 


OPTICAL  DENSITY/FT-  .565 
PERCfT  TRANSMISSICN/FT-  27.2 

OXYGBI  CCNCSSfrRATICN  IN  UPPER  ZCNE 
OXYGEN  PERCNT  3Y  VOLUME"  14.57 


CAS 


C02 

CO 

SCL 

HQ* 

HP 

NO 2 

SQ2 


CCNCafERATiaj 

(P«) 

.2469751E+05 

.7742417E+03 

.4374796E+03 

.1061680E+02 

.7624775E+02 

.5521900E+00 

Jl 

b . 

0. 

0. 

0 


LEVELS  EXCEEDED 
(0*0  1-YES) 
LVL1  LVL2  LVL3 
1 1 0 

1 1 0 

1 1 S 

110 
111 
0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 


DISTINCT  FIRES  AT  START  OF  FLAME  SPREAD  CALCDLAUCNS 


FIRE  NO 

ZCNE 

DIST-FHC 

BASE  FROM  FLOOR (FT) 

FLAME  HEIGHT  (FT) 

FIRE  BASE  AREAISQ  FT) 

BASE 

RADIUS. FLAME 

VOL  (FT) 

1 

LHR 

0.00 

3.25 

1.50 

.60 

2 

LHR 

0.00 

2.92 

1.75 

.35 

3 

LHR 

1.17 

4.24 

3.50 

.78 

4 

LHR 

1.18 

4.27 

2.75 

.79 

5 

IMR 

1.18 

2.57 

1.25 

.42 

6 

LHR 

1.18 

1.71 

.25 

.25 

ELEMENT  STATE  SUmARY  - 

OCNDITICNS  CN  ALL  SURFACES  AT  END  OF  FLAME  SPREAD  CALCDLATICNS 

1 2 

3 

4 ! 

5 

6 7 8 9 

10 

u 

12  13 

14 

15 

16 

17  18  19  20 

21 

22 

23 

24 

25 

26 

27 

SMOLDERING 

0 0 

0 

0 

0 

0 0 0 0 

0 

0 

0 0 

0 

0 

0 

0 0 0 0 

0 

0 

0 

0 

0 

0 

0 

FLAMING 

13  15 

8 

0 

0 

0 0 0 0 

0 

0 

0 0 

f 

0 

0 

0 0 0 0 

0 

0 

0 

21 

0 

0 

0 

CHARRED 

0 0 

0 

0 

0 

10  0 0 

0 

0 

0 0 

0 

0 

0 

0 0 0 0 

0 

0 

0 

0 

0 

0 

0 

FLAMING  AND  SMOLDERING  AREAS  BY  MATERIAL  TYPE 

MATERIAL  NO  1 

AREA  AFLAME  3.25 

AREA  SMLDRG  0.00 


2 3 4 

5.75  0.00  0.00 

0.00  6.00  0.00 


5 6 7 

0.00  0.00  5.25 

0.00  0.00  0.00 


C-27 


IBIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
PROM  OOP  I P URN  I SHED  TO  liDC  


TIME-  60.300  SEC  AFTER  IQJITICN 


DISTRIBUTION  of  elemental  states  at  end  of  flame  spread  calculations 


INTEGERS  CORRESPOND  TO  STATES  OF  INDIVIDUAL  ELQt- 

1- AFBIINT  state 

2- SMOLDERING  STATE 

3- AFLAME 

4- QJARRED 

5- HEATING, NOT  IN  CONTACT  WITH  FLAME 

6- HEATING, IN  CONTACT  WITH  FLAME 
7<M0LDERING,  COOLING 

AFT 


r 


THIS 
FROM 


furbished  TO  DDC 


TIME-  60.000  SEC  AFTER  IOJITICN 


FOR  SEAT  GROUPS J-  1-  4 

J-  5-  7 

J-  8-11 

J-12 

J-13-18 

J-19-21 

J-22 


CCSHIOJ,  BOTTOM 
BAOCREST,U*  REAR 
BACKREST, UPR  REAR 
BACKREST, TCP 
BJOCREST,  FRONT 
CUSHION, TCP 
CDSHICN, FRONT 


SEAT  GROUP  NO  4 


22 

ill  mini 

21 

nmmn 

20 

mm  mi 

19 

nmmn 

18 

mmmi 

17 

nmmn 

16 

minim 

15 

nmmn 

14 

nmmn 

12 

nmmn 

12 

nmmn 

11 

nmmn 

10 

nmmn 

9 

nmmn 

8 

7 

nmmn 

nmmn 

6 

5 

nmmn 

mmmi 

4 

nmmn 

3 

2 

1 

1111111333 

nmnin 

nmmn 

Variable  and  Unit 


Variable  Name 


Record  3 


Time  (seconds) 

Concentration  of  the  first 
toxic  gas  (ppm) 
Concentration  of  the  second 
toxic  gas  (ppm) 


TIME 

C0NC(1) 

CONC ( 2 ) 


Concentration  of  the  last 
toxic  gas  (ppm) 


CONC(NTXG) 


C.5  PROGRAM  STATISTICS 

DACFIR2  is  written  in  Control  Data  Corporation  (CDC) 

FORTRAN  EXTENDED  but  except  for  a very  few  statements  conforms 
to  ANSI  standard  FORTRAN  IV.  The  following  descriptive 
statistics  pertain  to  use  of  the  program  on  a CDC  6600  computer 
system  using  the  NOS/BE  operating  system. 

Core  Storage  - 130,000  words  (octal) 

Number  of  Source  Statements  - 3129 
Compile  Time  - 21  CPU  seconds 

Execution  Time*  - 248  CPU  seconds  for  175  sec  simulated 

time  with  an  integration  step  of  0.01  sec. 

* 

Execution  time  is  highly  dependent  on  the  integration  step 
size  used. 

C.6  PROGRAM  AVAILABILITY 

Copies  of  the  program  code  and  sample  input  data  may  be 
obtained  by  contacting 

Mr.  Charles  D.  MacArthur 
University  of  Dayton  Research  Institute 
300  College  Park  Avenue 
Dayton,  Ohio  45469 
Telephone  (513)  229-3921 
or 

Mr.  Charles  C.  Troha  ARD  520 
Department  of  Transportation 
Federal  Aviation  Administration 
Trans  Point  Building 
2100  Second  Street,  S.W.,  Room  1400 
Washington,  D.C.  20591 
Telephone  (202)  426-8416. 


C-30 


APPENDIX  D 

DERIVATION  OF  SEVERAL  RELATIONSHIPS  PRESENTED  IN  SECTION  2 

This  appendix  contains  the  derivation  of  several  relationships 
presented  in  Section  2. 

D.l  EFFECTIVE  VIEW  FACTORS  FOR  UPPER  ZONE  GAS  RADIATION 

Equations  (2-2a)  and  (2-2b)  contain  the  terms  F^  and 
the  effective  view  factors  for  radiation  from  the  upper  zone  gas, 
at  temperature  Tu,  to  the  solid  surfaces  in  contact  with  the  gas 
and  to  the  surfaces  in  contact  with  the  cooler  lower  zone  gas. 

To  compute  these  terms,  the  radiant  exchanges  among  the  major 
cabin  lining  surfaces  and  the  upper  zone  gas  were  evaluated 
separately.  Exchanges  between  the  ceiling  and  floor,  gas  and 
floor,  upper  walls  and  floor,  ceiling  and  lower  walls,  gas  and 
lower  walls,  and  upper  walls  and  lower  walls  were  considered. 

The  magnitude  of  each  exchange  was  determined  by  computing 
the  view  factor  in  each  exchange  for  the  range  of  upper  zone 
thicknesses  from  zero  to  the  full  cabin  height.  It  was  found 
that  the  gas  to  surface  exchanges  were  always  much  larger  than 
the  surface  to  surface  exchanges,  and  so  only  the  contributions 
from  the  gas-surface  exchanges  were  used  in  computing  Fu  and  F^. 

To  speed  computing,  linear  functions  were  fit  to  the  more  com- 
plicated view  factor  expressions  for  each  significant  exchange. 

The  resulting  expressions  for  the  effective  view  factors  are 

(a)  for  the  rectangular  cabin  cross-section 

f1  » 0.7587  + 0.02754L 

Fu  - 1.0 

and  (b)  for  the  circular  cabin  cross  section 

! 0.4782  + 0.0784L  for  L < 2r  - C„ 

0.1965  + 0.1457L  for  L > 2r  - CH 
Fu  - 1.0 
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where  L is  the  upper  zone  layer  thickness,  r is  the  cabin 
cross-section  radius,  and  CH  is  the  floor  to  ceiling  height. 

D.2  FLAME  ABSORPTION  COEFFICIENT 

Equation  (2-3)  relates  the  flame  absorption  coefficient,  k^, 
to  the  smoke  generation  rate  of  the  burning  material,  p" , and  the 
flame  height,  h^,  for  the  purpose  of  computing  the  flame 
emittance  and  thus  the  radiation  feedback  to  the  fuel.  This 
relationship  should  be  regarded  as  very  approximate  due  to  the 
overly  simplified  assumptions  involved.  The  motivation  for 
developing  Eq.  (2-3)  is  that  since  smoke  data  is  available  the 
relative  radiative  output  of  fires  on  the  different  materials 
can  be  judged. 

The  derivation  of  Eq.  (2-3)  assumes  that  the  flame  is  a right 
circular  cylinder  of  height  h^  and  base  radius  y . The  volume 
of  the  flame  is  then 

(D— 1 ) Vf  - *yo2hf 

To  compute  the  concentration  of  smoke  in  this  volume,  the  average 
residence  time  of  a smoke  particle  in  the  flame  must  be  known. 

If  it  is  assumed  that  all  smoke  particles  are  created  at  the 
flame  base  at  the  rate  pn  per  unit  area  (which  in  reality  is 
not  true)  the  number  of  particles  crossing  any  plane  parallel 
to  the  base  is 

{ D—  2 ) P - Ag>"  » fry  2 jb " - TW  2vP„ 

o o s 

where  v is  the  gas  flow  velocity  and  p is  the  particle  density 

3 f 11 

at  the  position  of  the  plane.  Thomas,  et.  al 1 J give  an 
estimate  for  the  gas  flow  velocity  at  a height  z in  a buoyant 
diffusion  flame 

(D-3)  v - 0.36  (2z  AT  g)1/2 

T 

where  AT/T  is  the  ratio  of  the  temperature  rise  in  the  flame  to 
the  ambient  temperature  and  g is  gravity.  The  flame  is  assumed 

^ Thomas,  P.H.,  R.  Baldwin,  and  A.J.M.  Heselden,  "Buoyant  Diffusion 
Flames:  Some  Measurements  of  Air  Entrainment,  Heat  Transfer,  and 
Flame  Mergin,"  Tenth  Symposium  (International)  on  Combustion,  The 
Combustion  Institute,  1965. 
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to  be  at  a constant  absolute  temperature,  Tf  = T + AT,  and  the  gas 
at  the  flame  base  is  assumed  to  accelerate  from  zero  velocity 
so  that  the  number  of  particles  in  the  flame  volume  N is 


(D-4) 


/ p dV  » f 
Vf  ° 


Substituting  (D-3)  into  (D-4)  and  assuming  that 

AT/T  = 2000  °R/500  °R  = 4,  the  integral  can  be  evaluated 


giving 


(D-5) 


2£>"(hf/g)1/2A 


Now  making  the  further  simplifying  assumption  that  the  N 
particles  are  uniformly  distributed  in  the  flame  the  mean 
particle  concentration  is 


(D-6) 


C = JL 

Vf 


(h*g) 


The  absorption  coefficient  is  employed  in  the  Bouguer -Lambert 
law  (uniform  grey  gas)  as 

(D-7)  i * iQexp(-kfs) 

where  i is  the  intensity  at  a distance  s along  a straight  ray 
path  away  from  the  point  where  the  intensity  is  iQ.  The 
definition  of  a "particle"  of  smoke  is  based  on  the  Bouguer- 
Lambert  relationship  but  in  a different  form 


(D-8) 


"o  10 


-0 . 04575CS 


By  equating  (D-7)  and  (D-8)  C is  related  to  k^  to  obtain 
Equation  (2-3) 

kf  - (0. 105)C  « 0.21  £"/(hfg)1/2 


i 


